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ABSTRACT 


This thesis considers the problem of determining 
optimal product mix and machine loading for single and multi- 
stage production systems. It is assumed that the Group 
Technology (GT) concepts have been utilized to classify the 
jobs under consideration for manufacture on the production 
systems into groups. Each job has a lot size comprising of 
identical units and each production stage has onlA/’ one 
machine. In case of multi-stage production systems, it is 
further assumed that the job operations do not have any 
technological ordering. Corresponding to each stage, the 
job available for processing either have prespecified 
production times or the production times are to be determined 
optimally considering machining speeds as decision variables. 
The optimization is to be performed considering maximization 
of production rate or maximization of total profit. 

For each one of the above stated cases, an optimizing 
algorithm which uses the concepts of branch and boiind proce- 
dure is developed. 

In order to save on the computational effort involved 
in the use of branch and bound procedure, for each case a 
heuristic method has been developed. The heuristic method 
helps in reducing the size of the problem to be handled by 
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the branch and boimd procedure. Illustrative examples are 
given to explain the various steps of the optimizing and 
heuristic algorithms. 

The optimizing and heuristic algorithms have been 
computerized and implemented on DEC--1090 Computer System, 
Based on the limited experience of solving the numerical 
examples? it was observed that the heuristic algorithms have 
computational superiority over the optimizing algorithms. 



CHAPTER I 


INTRODUCTION 


1.1 GROUP TECHNOLOGY: 

Group Technology is a very progressive method of 
organising production and especially it is becoming very 
popular in those industries which are engaged in medium and 
small batch productions. Group Technology concept identifies 
and makes use of sameness or identical nature of operation 
processes and parts in design and manufacture. Similar parts • 
are grouped together into one part family based on similarity 
of their design and dimensions, geometrical shapes and technolo- 
gical requirements, and according to the needs to process 
designated part families, machine group or cells are formed. 

1.2 REVIEW OF PREVIOUS WORK; 

Analytical methods suggested by Burbyidge [l] take into 
consideration those informations which are contained in route 
cards. However, he does not consider design and geometrical 
similarit 3 r. Batra and Rajgopalan [2] also use informations 
given in route card for each job and using graph theoretic 
approach they form machine cells. McAuley [3] uses similarity 
coefficient for every pair of machines to form group of machines. 
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El-Essa'wy[4] does component flov/ analysis to form machine and 
component grouping. 

Many researchers have proposed models to schedule 
products on machines available in machine cells. Petrov [5] 
has developed four interrelated scheduling models for different 
type of route sequences and com.ponent flows. This approach does 
not guarantee optimality mathematically. Ham and Hitomi [6] have 
suggested a technique to schedule various grouped products on 
multi-stage production system. 

In addition to the problem of group formation and group 
scheduling, machine loading and product mix decisions represent 
major problem areas for group production planning and schodul- 
ing. Machine loading problem determines the total number of 
units of various jobs to be processed within limited capacity 
of various production facilities available, while product mix 
decision determines the kinds of jots to be processed within 
limited production time. Hitomi and Ham^ [7,8] have considered 
the problem from the view point of Group Technology for a single 
machine production system. They have presented models for' two 
distinct situations. In the first case, the unit production 

i 

time is prespecified [7], while in the other case, it is 
dependent on the machining speed which is treated as a decision 
variable [8]. They maximize total number of units produced of 
various products and in process do not account for weightages 



or importance of various products to be produced. Moreover, the 
experience of the author in using Hitomi and ham' s algorithm on 
a number of problem suggests that their algorithm does not 
yield optimal solution necessarily as claimed b 3 r them. 

1.3 SCOPE OF PPuESEKT STUDY? 

Basically, this thesis is concerned v/ith the machine 
loading and product-mix decision problem considered hy Hitomi 
and Ham [7,8] with added features to make the problem more 
realistic. In real life, single stage production environment is 
not very common. Further, besides maximizing the total amount 
of production, the management ma^r desire to evolve a product mix 
and machine loading stratog^^ which maximizes the profit. There- 
fore, in the present work, an attempt is made to develop models 
and solution methodologies for multi-stage situation with 
maximization of total production or maximization of profit as 
the criterion of optimization, considering the unit production 
time as prespecified or dependent on the machining speed ( a - 
decision variable ). In general, the solution methodolog 3 ^ 
exploits the branch and bound concepts for its devoloj3ment. 

Since, the use of branch and bound approach involves conside- . 
rable effort, heuristic based computationall 3 ^ efficient solution j 
methodologies which initially reduce the size of the problem 

I 

and then find the solution applying branch and bound concept | 
to the largely reduced problem, are also presented. i 
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Further^ the algorithm for single machine loading 
problem and product mix decision of Hitomi and Ham [7,8] has 
been modified. The modified algorithm ensaros optimal solution 
to the problem. 

1.4 ORGANISATION OF THESIS; 

The cases of single stage production system have been 
dealt in Chapter II, In this chapter besides the modified 
version of algorithm of Hitomi and Ham [7,8], optimal and 
heuristic solution procedures for maximization of profit are 
also presented. Various variations of the multi-stage produc- 
tion system have been considered in Chapter III. Numerical 
examples are presented to illustrate the solution methodologies 
proposed for each of the variations of the single and multi- 
stage production systems considered. In Chapter IV, conclusions 
based on the present study alongwith suggestions for further 
work are presented. 



CHAPTER II 


SINGLE STAGE PfiOBLEJ'I 

2.1 STATEf'IEOT OF THE PROBLEH: 

Consider a single stage production system with limited 

available production timeoN types of jobs are available for 

production. Based on part-family concept of group technology, 

these jobs are classified into N groups. Let groups 

e 

Gj^ (i = 1,2,.,., Ng) contains Ih jobs. 

The objective is to determine an optimal machine loading 
and product mix such that within the available production time, 
the quantities produced of each job do not exceed the lot size 
limit. The optimization is carried out to maximize either the 
production rate or the total profit. 

Corresponding to each group, there is a group production ; 

t 

time which is comprised of the group set-up time and the sum ■ 
of the job production times of the jobs in this group. The job 
production time includes the job set-up time and the time 
required to produce the various units of that job. The unit : 

I 

production time required to produce single unit of any job ma,y 
be prespecified or dependent upon the cutting speed (a decision I 
variable) . 
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2.2 ASSUMPTIONS? 

The various assumptions are; 

i) The single stage production system contains onl3r one 
machine . 

ii) Only one unit of the job can be processed on the machine 
at a time. 

iii) A limited amount of time is available for production on 
the machine. This is the actual productive time available 
on the machine and as such do not include time lost due to 
breakdowns, etc, 

iv) Each unit of the job must be processed to completion. 

v) All the units of a job are produced together. Similarly;, 
all the Jobs belonging to same group are produced together. 

2.3 NOMENCLATURE? 

Following notations are used for the development of the 
mathematical models? 

a. . preparation time for job J. . (min/pc) 

-i- J 1. J 

bfj tool replacement time for job (min/pc) 

B. . profit per unit of job J. . (Rs./pc) 

J 1 J 

C. . 1-minute tool life machining speed for job J. . (m/min) 

-L J 3. J I 

^u total capacity of the machine utilized (min. ) 
d available production time (min) 

E. . high-efficiency machining speed range for job J. . (m/min) 
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N.) 


"id 


10 

^id 


N 

N. 


N, 


P 


id 


(c) 


i-th group or part family 
group index ( = ls,2,...s N^) 
job index (= 1,25.. 
j-th job in i-th group 
lot size of J. . (pcs) 

_L 

slope constant of Taylor tool life equation 

total number of jobs available for production 

total number of part-families or groups 

total number of jobs in group 

unit production time for J. . (min/pc) 

10 


p. minimum-production -cost unit production time for job 


10 


(t) 


(min/pc) 


p. . ' maximum-production-rate unit production time for job 


10 


J. . (min/pc) 

J 


p. .(v. .) cutting speed dependent unit production time (min/pc) 

P. . job production time for J. . (min) 

1 J 1 J 

P- ^ maximum™.production_ra.te job production time for J. . (min) 

1 J J_ 

q^^(v. .) speed dependent unit production cost (Rs/pc) 


"id 


10 

T. . 
10 


U 


^id 


job set-up time for J. . (min) 

X J 

Group set-up time for group G^ (min) 
actual machining time for J. . (min/pc) 

1 J 

tool life for J. . (min/edge) 

i J 

set of all parts available for production 

l»2j..«j N } 


, N. ; 1 


{ d “ lj2,... 

machining speed for J . - (m/min) 

X J 
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V . . 

ID 

^iD 

^iD 

w. . 

ID 

^iD 


Y 

y' 

Z 

z 

z 

z. 

a 


Cc) 

(c) 

(t) 


minimum-p reduction -cost machining speed for J. - 

J- J 

high efficiency machining speed for J. - (m/min) 

i- J 

maximum-production rate machining speed for J. 


ID 


max 


^iD 

Y . . 
ID 

6 

Aid 

0 


quantity of to be produced (pcs) 

0-1 type variable for J. - 

ID 

0-1 type variable for group 

total production cost (Rs. ) 

total variable production cost (Rs.) 

value of objection- function 

lower bound on objective function value 

upper bound on objective function value 

maximum value of objective function 

direct labour cost and overhead (Rs./min) 

mochining overhead for J . - (Rs./min) 

ID 

tool cost for J. . (Rs./edge) 

_L J 

slack time (min) 
machining constant for J . . 

-L J 

null set 

lagrange multiplier 
shows optimal value 


Following additional notations will be used in 


(m/min) 

(ra/min) 


solution procedure: 
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D consumed production time (min) 

n all nodes belonging to N1 

N1 set of all descendants 

R set of remaining jobs at any node 

s time still available for proiuction (min) 

S set of selected jobs at any node 

shows candidate for acceptance or rejection 
[a] is a Guassian notation imploring a greatest integer less 
than or equal to A. 

2.4 MATHMATICAL MODELS AND OPTIMAL SOLUTION PROCEDURES J 

Mathematical models for the various cases of the single 
stage problem stated in Section 2.1 are presented alongwith the 
optimal solution methodologies. 

2.4.1 Ca se 1: Maximization of Total Units _Pugdug _s5l ' 

In this case wo maximize the total amount of production 
irrespective of the relative values associated with the various 
jobs to be processed. Two variations of the basic problem are 
considered. In the first variation we assume that the unit 
production time is prespecified while in the other case v/e 
assume that the production times are dependent upon the machining 
speeds. The machining speeds are also to be determined opti- 
mally considering the workpiece - cutting tool-machine tool 
considerations. 



2.4. 1.1 Model 1: Unit Production Time is Pre specified? 
Constraints ? 

Various constraints on the problem are listed below, 
i) Constraint on Machine Utilization? 


This constraint ensures that total time consiimed for 
production does not exceed the available capacity of the machine, 
Mathematically, this is represented ass 




1 


E hd "id^ ^ 

= 1 J-E 


( 2 . 1 ) 


where , 


^ij production time for producing w^^^ units of 


job It is expressed as. 


hd 


s. . + w. . p. . 
10 10 10 


( 2 . 2 ) 


_ fl if jot) Jii is selected for productionwp 

ij “ ‘0 if job J. . is not produced *2^ 

1 0 


and, 





if 


1 

Z X. . > 1 
0=1 “ 


otherwise 


(2.4) 


ii) Constraint on Units Produced? 

The quantity w. . to be produced of job J. . should not 
exceed the lot size specified for it. Thus, 
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O bjective runction ; 

Since the objective is to maximize the total number of 


the units produced, we maximize a function presented asi 

N N. 
g 1 

Z = Z Z w. . X. . 
i=l j=l 


( 2 . 6 ) 


Solution procedure s 

Let Q represent the time needed to produce all the jobs 
to their lot sizes. Then, 

N, N, 

(2.7) 




^0 ij -10' 

If Q < djthen all the jobs will be produced equal to their lot 
sizes and the optimal solution will boi 


where. 


■¥r 

^id 

= 1 

¥ J. . e U 

10 

(2.8) 

w? . 

= 1. . 

¥ J. . e U 

(2.9) 

Z* 

= z 

max 


(2.10) 

z 

max 

N N. 

g 1 

= Z Z 

i=l j=l 

^id 

(2.11) 


If Q > d, then some of the jobs are either not accepted for 
production at all or are not produced to their lot sizes. In 
such a situation, obviously, 


Z* < 


max 
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A branch and bound procedure is developed to optimally 
ascertain the '■machine loading and product mix for such a 
situation. 

In the branch and bound procedure, at any stage we branch 

rsj 

a node n for which the upper bound value Z, is maximimn among 

all descendents. Z is calculated as the sum of the lower 

S.(l-X. )+s- . 

bound value Z and maximTom of + p. . }] among 

ij 

all remaining jobs at that node? where s shows the remaining 
time available for production. Initialises, Z = 0. If at this 
branching job Jrv is selected for production, the lower bound 

X J 

value Z, is incremented by wr-r. 

Thus, the branching procedure finds the lower bound and 
the upper bound values for each of the descendents. The node 
with maximum upper bound value is identified. For this node, 
if the upper bound value is equal to the lower bound value, 
then branching procedure is terminated and this node carrying 
the maximum upper bound value gives the optimal solution. 

Optimizi ng algori thm s 

maeaAssr tti 'a n - rn ; :; mi rTiTrmir r rihiii-niiir imm i n iii|im..i u i[aii i L.iiiimiiiimwLii 

The following is the step bsr step procedure for deter- 
mining the optimal machine loading and product mix. 

Step I s Calculate Q, from (2.7). If Q <_ d, then go to Step 2J 
■otherwise, go to Step 3. 

Step 2 ; Manufacture all the jobs to their lot sizes. The optimal 
solution is given by the following relationships. 
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Terminate the procedure. 

Ste p 3 i Follow Branch and Bound Procedure. ; 

Select an initial node; set n = 1, N1 = {n}k For this 

node n, set R = Uj, S = 0s, s = d 5 Z»0 and = 1 ¥ i. 

Find a nob J-r-r from the nobs J. . e R for which 
^11 '-'10 

(Sr + ST-r)/l-r-r + p-r-T is minimum. Find 
Z = [d/{(Sj + + Plo^^ 

Step 4 ; Find a node from the nodes n e N1 for which Z is maximum. 

Denote this node by n and J^r-t in n by Jc-y. For the node 
n, if Z = Z, go to Step 85 otherwise, go to Step 5. 

Step 5: Find D = D + S'xXr + src' + wc*r p-rc', where wcp' = max w,. 

— 11 10 10 10 10 10 

w. . is varied from zero to lr~ till D <_ d. If wrv = 0, 

go to Step 65 otherwise, go to Step 7 . 

Step 6 ; Remove the oo'b from the set R of the node n. Set 

R = R~{j~~}. If R = 0, set Z = Zj otherwise, find 

s'?X. + s. . 

min ( - 1 - p ) and denote this oo'b hy O’jy 

Jp^eR -^id 

Set Z = Z + [s/KS-Xj + + Pj-B. 

Return to Step 4. 

, jfSJ 

step 7 i Branch the node n into two nodess 

First nodes Set n = n+1 and N1 = N1 + {n} - {n}. Then 
for the newly created node n, R = R - { J'X'y} , S = S, 

X J 

D = D, s = s, and 2 = Z. If R = 0, set Z = Z; otherwise. 
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find min {(S. X. + +P--;} and denote the 

J. -eR ^ ^ 13 13 

13 


corresponding job by Set ^ = Z + [s/{^jXj + 


+ PI3 }]. 


Second nodes Set n = n+1 and N1 = N1 + {n}. Then for 

the neyly created node n, R = R, S = S + {3''r~} ,0 = 0, 

1 J 

s = d - Dj, Xo' = 0 and Z = Z + w'y'v. If R = 0, set 

1 "**" "*** 1 o 

Z = Z| otherwise, find min {(S.X. + s. .)/l. .+p. .} . 

T pR ^ 

ij 

Denote the corresponding job by and set 

Z = 2 - [s/lCSjXj . SJ3)/1J3 . Pj,}]. 

Return to Step 4. 


St ep__8 s Optimal solution will be contained in node n, Set, 


,* 


,* 


3 = S, Z = Z and stop. 


The value of w^ required in Step 5 is determined 

J. J 

using the following relationship. 




ci^D ^ SrXr 

pr-:r ^ 
10 




ij.] 


If wro" > Irn", we set wrc = Irc^. 

13 13 13 13 


Nu merical exampl e! 


The example given is based on the data of a single 
stage problem considered by Hitomi and Ham [7]. 

Consider a single stage problem where ten jobs are 
available for processing. These jobs are divided into four 
groups. The relevant data for the problem are given in Table 2. 



15 



Table 

2.1s Basic 

: Data for 

Numerical 

Example 1. 

GROUP NO. 

PART 

LOT SIZE 

GROUP SETUP 

JOB SETUP 

UNIT PRODUCTION 


NO. 

pcs 

TIME (min) 

TIME 

(min/ lot) 

TIME (min/pc) 


1 

30 


19 

6 

3 

2 

40 

40 

8 

2 


1 

20 


10 

17 

2 

2 

50 

35 

9 

9 


3 

30 


15 

7 


1 

20 


5 

12 

3 

2 

60 

20 

13 

6 


1 

20 


6 

16 

4 

2 

10 

45 

10 

15 


3 

40 


20 

13 


Available Time d = 40 hours. 


Solution s Using the various steps of the algorithms, we 
obtain, 


Is Q = 3105 mins. 

Since Q > d, go to Step 3. 

Step 3 i n = 1, N1 = {n} = {l}. For node 1, set 

R = ^21' *^22’ '^23* '^ 31 * '^32’ *^41* ‘^42’ '^43^" 

S = 0 , s = 2400,00 mins. , D = 0, Z = 0 and = 1 ¥ ±. 

Calculate {(S. s. -)/l. . + p.^} for J. . e R. The 
minimum value of 6,55 min/pc corresponds to ;iob J^2* 
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Step 4 ; 




Thus, Z = 366, Proceed to Step 4. 

Since for node 1, Z ^ Z, go to Step 5. 

For node 1, D = 393 mins, and Therefore, 

go to Step 7. 

Branch out from node 1, two nodes, 2 and 3- 
At node 2: Set n = 2, N1 = {l} + {2} - {l} = {2}. 

For the newly created node 2, set R = R - '^‘^32^ 

= Ji2» '^21’ "^22' ‘^23’ ‘^31" "^41’ "^42’ ^^43^' 

s = 2400 mins,, D = 0 and S = 0. Since R ^ 0, 

calculate (S-X- + s. .)/l. ^ + p. . for e R. 

Identify the job for which the value of this expres- 

sion is minimum. For the present problem the minimum 

value of the expression corresponds to job J22* 

Similarly, at node 3: Set n = 3, N1 = {2} + {3} = {2,3 

For the newly created node 3, set 

R = R = J41sJ42»J43} » 

D = 393 mins., s = 2007 mins, S = ^‘^32^’ ^ ~ 

= 0. Since R / 0, we find the minimum of 

(S.X. + s. .)/l. - + P-^. Its value is 6.8 min/pc 
corresponding to job J2.2* find that Z = 355. 

Return to Step 4. 


From Step 4, we find that max Z occurs at node 3. 

neNl 

Further branching takes place at node 3. The steps 4 to 7 
of the algorithm are repeated till the branching procedure is 
fathomed. A portion of the tree-diagram obtained using the 
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A portion of the Tree-Diagram obtained for Numerical 
Example 1 using Optimizing Algorithm. 





17 


branching scheme for this example is depicted in Fig. 2.1. 
The optimal obtained is given in Table 2.2. 


Table 2.2: Optimal Solution for Numerical 
Example 1. 


Index of Accepted 

Quantity?- to be 

produced 

Group 

Part 



1 

1 

30 


2 

10 


o 

2 

50 



3 

30 


~z 

1 

20 


D 

2 

So 


4 

2 

3 

7 

40 


Total 

number of units 

produced '= 

247 

Total 

time utilized 


2390 mins. 


2. 4. 1.2 Model 2: Machining Speeds are Decision Variables j 

Before structuring a mathematical model of the problem, 
we develop mathematical expressions for job production time and 
cost of production. 

Job Production Time : 

The unit production time p. .(v. .) for job J . . is a 

function of machining speed v. . for that job and is given by, 

^ j 
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* *13 * *13 ^ 


or. 


p. .(v. .) = a. . 


(i/n . -1) 

:x . . b . . V. . 

10 10 

17H77 

"ii 


iO 


¥ J- . e U 

■1- tJ 


(2,12) 


where a. . is the preparation time, b. . is the tool replace- 

ment time, t. . is actual machining time, T. . is the tool 
■^O -i- J 

life,7\^^ is a machining constant, n^. ^ and C^. ^ are parameters 


11 


11 


for the Taylor tool life equation for oo'b 

1 J 

Let the unit production time given by(2»12) be 

(t) 


minim'ura corresponding to a point v. . 

X J 


This machining 

speed referred to as the 'Maximum- Production- Rate Speed' or 
* Minimum-production- Time Speedy is obtained by setting first 
derivative of p. .(v. .) with respect to v . . equal to zero. 
This yields. 


V. . 
10 


(t) 




n. . 

11 


¥ J. . e U 
11 


(2.13: 


If w. . units of J. . are produced at a speed v. 

J- ^ X 


the oot* production time, P ,>^9 is given by, 

J 

h3 = ^13 


( 2 • 14 


Production Costs 


The per unit production cost, q- ^(v. .), is also a 

<J X vJ 

function of machining speed v, . and is expressed ass 

X J 
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qy(v.j) = a a. J (a + h- (<xb.j ■. Y. .) 5;:i^ 


or, 


pi • ' 

q, .(v, .) = aa, . + (a+p,.) 


10 


10 


10 V 


10 


^io ^ 


(1/n ~1) 


¥ J. . e U 

J> J 


(2.15) 


where a is the direct labour cost, is the machining 

overhead and y . . Is tool cost, respectively , 

10 

f C ^ 

Let V. represent the machining speed correspond- 

-L J 

ing to the minimum cost point. The value of this point is 
obtained by setting the first derivative of (2.15) with 


respect to v. . equal to zero. The cutting speed, v. . 


(c) 


10 


10 


IS 


called the ’ Minimum-Production-Cost Machining Speed ' and 
is expressed as follows s 

(c) _ 


"^iO 


(a + 

c { 


n. . 

} 


(l/n.^-l)(aby 


¥ J. . £ U (2.16) 

i J 

Constraint t 

The various constraints for this model -are the same 
as that for Model 1. The constraint set is given by (2.1) 
and (2.5). 
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Function s 


In this model we have two different objectives to 
consider. The first objective is to maximize the total number 


of units produced \\Ahile the second objective is to minimize 


the total cost of production. These objectives are referred 
as the primary and secondarA;" objectives, respectively. The 
expression for the total number of units produced is the same 
as in Model 1 and is given bj?' (2,6). The expression for the 
cost of production Y,is given ass 


Y 


N 

g g 

I {a S.X, + E 
i=l ^ j=l 


(a6 . . 


“id ‘iid<nd» ^id^ <2- 


Optimization is carried out using the solution proce- 
dure which initially maximizes the primary objective function 
and then it minimizes the secondary objective function. 


A nalysis ; 


Initially for all jobs J. . e U, we calculate the 
minimum-production-time machining speed, v. from (2.13) 

and the unit production time is then calculated from (2.12). 

The following expression is used to calculate the job production 
time for each job considering that each unit in the lot of that 


job is produced at minimum unit production time; 

P. = s. . + 1. . p. . (v. 

ij 10 10 ^10 ^ 10 ^ 


17 


( 2 . 18 ) 
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Let q"^ represent the total time needed to produce • 
all the .jobs to their lot sizes. Note that each unit of the 
lot is produced at the minimum-production-time machining 
speed. Mathematically, 

t % / (t)\ f \ 

or = s (s. + I p. (2.19) 

i=l ^ j=l 

Depending upon the value of q"^, there are three possible 
situations. These are; Q**" = d, q"*" > d and q”*" < d. These 
three cases are discussed separately. 

+* i" 

■i) Q = d ; If the total time Q , needed tc produce all the 
units is equal to the available capacity d, of the machine, 
then all the jobs to their lot sizes are produced at the 
minimum-production-time machining speeds. Therefore, the 
optimal product mix and the cutting parameters will be given 
by the follov/ing relationships; 


= V 

^ij ’ 

e U (2.20) 

where U is the set of all the jobs available for processing. 


Z* = El. 

J. .eU 
ij 

13 

= 1» V. . 

* 13 

w.*. = 1. . 
la 13 


^ J. ^ 
13 


i^) ^ > d ; When > d, obviously, all the units cannot be 

produced even when minimum- production-time machining speeds 
are used. Any one of the following situations may result. 


1. All the jobs are taken up for production. However, 

for one of the jobs all the units cannot be produced. 
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2. All the oohs are not taken up for production. The 

jobs taken up for production are either produced to their 
lot sizes or for one of the jobs, the desired let size 
cannot be produced. 

The optimal product mix and machining speeds cun 
be obtained using the following approach. 

First of all, v;e determine the optimal machine 

loading and product mix based on the minimum-production- 

time speed for each job using the solution methodology 

discussed in Section 2,4. 1.1 for Model 1. Let S* and w?. 

ij 

(J-. e S*) respectively represent the set of jobs accepted 

for production and the quantities of each job to bo produced. 

Then, the optim.al solution turns out to be the following; 

1 ¥ J. . e S* 

X* = { (2.21) 

0 Otherwise 
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If slack time 6 ~ 0 , i.e. if the machine is utilized 


to its full capacity, then all the selected are produced 

at the minimum-prc duct i on-time machining speeds. Obviously, 


v3. = V 

10 iO 


^ J, . s S 


(2.25) 


If 6 >0, then we have positive slack time which, can be 
advantageously utilized to minimize the total cost of 
production. It is obvious tliat the minimum-producti on-time 
machining speed will be higher than the minimum-producti on-cust 
machining speed and tht^ cost of production at the minimum- 
production-time machining speed will be higher than that at 
minimura-production-cost machining speed. So, it is reasonable 
to assume that the unit production cost is least at the miniraum- 
production-cost machining speed and it increases with an 
increase in the machining speed upto the miiiimum-production- 
timo machining speed. Ther''fure, the slock time can be utilized 
in lowering down the cost of production by decreasing the 
machining speed from to v. The objective is to 

determine v*., v. < vf. < v. 


, such that the total 
cost of production is minimized. The problem can now be for- 


mulated as a nonlinear programme given below: 

N, 


g 


N. 

1 


Minimize Y = I (aS.X. + E (as.. 

^ ^ j=l 


i=l 


+ w'?. q. . (v. .)) X*.] 
n ii 11 11 


( 2 . 26 ) 


subject tc, 
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N 


N, 




The inequality given by’ (2*27), represents the machine capacity 
constraint* The Equations (2.26)and (2.27) contain certain 
constant terras which can be ignored for the purpose of 


optimization. ’ 'The optimization problem; can be restated asy 

N N. 

g jL 

Minimize Y’ = 2 S w* q . . (v. .) X* ' , 


(2.28) 


subject to, 

N N. 

2 w* . p. . (v. .) X* . < d 

i_l i_i 11 11 11 11 " ^ 

JL J., ^ 


where , 


2 

i=l 


N. 


:Y - 2 (a S . + 2 a s . . x. . ) 

1 1 11 11 


(2.29) 


(2.30) 


and 


N. 


N. 

1 


do = d 


2" (S.X* + 2^ s, . x!. ) 

1 1 11 11 


i=l 


(2.31) 


For the problem represented by (2.28) and (2.29): 
a Lagrangion function of the following form is obtained) 


I* (v,. ^ , p ) = 


N N. 

g 1 


11 


i=l j=l 


2 w*. q. . (v. .) xT. 


10 ^10 '10 11 


H N. 


+ p { 2'" 2 wt. n. . (v. .) x?. + Ft 

i=l 0=1 






(2.32) 


where p and k represent the Lagrangian multiplier and slack 
variable, respectively. For optimal values of vf . , p' and k^ 

X J 

we set 9 l/ 3 v^j = 


0. It gives. 
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or, 




-3 


3v;. , 


■¥r 

ij 


Otv* . 

^ 1.1 


= 0 


a+P . . 

■( ~i )Aii 


_iJ. 


V*. 

11 


l/nT. 


11 '11 -^*11 n -'-11 

4i 


, (l/n,^»2) 


"ii 


) = 0 


From the above equation, we get, 


1 "Sf 


P • - - ( 

11 


V. 


^ 1 ) -Y ^ 

J n 


l/n. . 
11 


i-i 


-1 ^ ^ ’’id ( "4d' 


V*. l/n,. 


11 


a 


(2.33) 


^ a + ii'* + P 


n. 


¥ J., e S^ 

11 

(2.34) 


Similarly, setting. 


'd^ 


0 yields. 


N 




= d - S E . p. . (v*. ) X* . 

o 3^1 11 ^11 ' 11 ' 11 


(2.35) 


and ^L/gk = 0 gives. 


p*k*= 0 (2.36) 

If we substitute v? . = v. in (2.33)? we find that 

IJ 1^ /.N 

= 0. Similarly, if we substitute v?. = v. in (2.33) 

J- J ” a- J 

we obtain, = “=>. Hence, It is obvious that, 
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> 0 for v’r. e (v v (2.37) 

-1-0 i U «L o 

Since (i ' is nonnegatives it guarantees the minimization of total 
production cost. 


If n' 


V. 


13 


Oj then from (2.36), k’ ^ 0 and from (2.34), 

(c) 2x T . a* (2.38) 




¥ J. . e S 
10 


However, if > 0, then from (2.36) v/e knov/ that k"' = 0. 

The optimal values of ji* and vf . e E. . can now he determined 
from (2.33) and (2.34), respectively, by using hit and trial 
approach. E. . is called the high-efficiency speed range. 

J_ J 

Mathematically^ , 


E. . = [v. . 

10 10 


(c) 


(t) 


id 


] 


(2.39) 


iii) Q < d ? When Q < d, the optim.al solution turns out to 
be the followings 


R* = 0, S* = U, X* = 1 ¥ is 2* = 1. 

" 1 “<1 «« = ha ^ ha ^ s* (2.40) 

In this situation the slack time is expressed as, 

6 = d - (2.41) 


We can minimize the total cost of production by 
utilizing the available slack time. Once the available slack 
time is obtained, the optimal machining speeds can be determined 
by foraulating a nonlinear programme of the problem given by 
(2.26) and (2,27). The optimal machining speeds will be given 
by (2.34). 
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mizing Algorithms 


A step-by-step algorithm, for finding the optimal 
prciduct mix and machine loading when the machining speeds 
are also decision variables, is given below: 

Stop 1: For each .lob J. . e U find v. and using 

j 

( 2 . 13 ) and ( 2 . 12 ) , respectively. From (2.18) and 
( 2 . 19 ), find Q*". If q"^ £ d, go te' Step 2j otherwise, 
go to Step 3 . 

Step 2: The optimal solution is. 


0 = 

U 



X. . = 
10 

1 

¥ 

‘^id 

wf . = 
10 

do 

¥ 

'^id 

xt = 

1 

¥ 

i 


J.-eS'^ 

10 

If o'" = d, set v/. = V. ¥ J. . e S*^ and stop. 

For q"^ < d, go to Step 9. 

Step 3 i Follow branch and bound procedure. Select an initial 

node: sot n = 1, I'll = {n}. For this node, Set R = U, 

S = 0, D = 0, s = d, Z = 0 and Xj_ = 1 ¥ i. Find a 

job Jt-v from the jobs J. . s R, for which 
10 10 

(S:|- 4- s^-;^)/lj"^ + p:r-;T is minimum. 


Set Z 
Step 4. 


[ 4/ {(Sr + srv)/l^ + prvj] and go tc 
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Find a node n from the nodes n e Nl, for -which Z 


is maximum. Denote the job in n by For the 

1 J 1 J 


node n, if Z = Z , go to Step 85 otherwise* go to 


Step 5. 


Find D = D + Sr Xr + sr:p'+ wrr prr 


1 1 


la "la 


where 


'w'yy = max w. w. . is varied from zero to Irr till 
10 10 10 ^0 


D £ d. If wrv = 0, go to Step 6; otherwise, go to 
““ 10 


Step 7. 


Remove the O'^'b Jrv from the list R. Set R = R-{J^}. 

X -i- 


If R = 0, set Z = 2. Otherwise, find 


min {(S.X. + s. .)/!.. pv"<}, denote the corres- 

J £ R -X X XJ xj 

ponding Job) by and set Z = Z + [s/{ (S'^Xj+S'rrT)/!^’’ 

Xj XX Xtj xj 

+ pf^}]* 


Return to Step 4, 


Step 7 s Branch the node n into two nodoss 


At first nodes n = n+1, Nl = Nl + {n} - {n}. 


For node n, set R = R - {Jrr}, S=S, D=D, s=s 

Xo 


and Z = Z. If R = 0, set Z = Z. Otherwise, find 


min {(S.X. + s. .)/l. - + P -V } • Denote the corres- 
J eR 11 10 10 10 

iO 

ponding job) by J*. tir and set Z = Z + [s/ICS^Xj + 

XJ XX ij XJ 

+ p^V}]. Similarly at second node: n = n+1, 


iCL = Nl + {n} . For this node n, set R = R, 


S + {Jr'^ 


D = D, s = d - D, Xr = 0 and Z = Z + wm^y. If R = 0, 


set Z = Z, Otherwise, find 


d min { ( S . X . +s . J /l . • ^-p - . ^ ^ 


J . -eR 

XJ 


X 1 XJ'^ XJ 
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Denote the corresponding job by Jrr and set 

Z = Z * [s/{(s-x, + SJ3)/1J3 . 

Return to Step 4. 


U The optimal solution given by list n, is. 
„ * 


= S 


X. 


* 

'id 

* 


1 ^ J. . e S 

■0 otherwise 


X.. 


* 


0 if .1. xd . = 0 

; 1=1 11 
'• 1 otherwise 


Go to Step 9. 

Find from (2.16) and from (2.12), 

respectively. Calculate 


N. 


N. 


Q = {S^X^ (s.j + w. . p.^( >) Xy } 

If Q < d, set vp = V. ^ £ S* and stop. 

Otherwise, go to Step 10. 

Step 10 ? Select a job J. - e S*. Find initial value of |j. from 

J- ^ 


the following expression, 

P • ■ - (" ^ 

^11 1 

4 = — 


V. + V. 

1) Y . ( -iJ f.T 

n. . id 20. . 

11 11 


l/n. . 

) 


.1 H. ( -i_ _1) b ( “ 

"io 


mr 




- a 


Go to Step 11. 

f p) 

Step 11 s For the value of [j, obtained in Step 10 , find v . o: 

J- J 

and p. . from (2.34) and (2.12), respectively. Calculate, 
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Q = 


N 

g 

I 

i=l 


{S.X. 

^ 1 X 


* ■"« <d * 


If Q ji^ d, go to Step 12? otherv-zise, stop. The set 
of machining speeds obtained currently is optimal. 
Step 12 ; Modify the value of ^ appropriately depending on 

whether Q < d or Q > d. If Q < d, reduce the value 
of p and if Q > d increase the value of p . 

Go tc Step J.I.. 

A systematic procedure for the selection of the value 
of p is presented in the form of a block-diagram in Fig. 2,2. 


riiap '. erical Example 2; 


For better understanding of the solution procedure, 
consider the single stage problem considered by Kitomi and Ham [8]. 
In all, there are ten jobs available for processing. These jobs 
arc divided into 4 groups. The basic data are given in Table 2.3- 


Sol\itions 


Initially, the values of v. and p. are found 

_L ^ O 


for 'all jobs J. - e U. These are shoT.m in Table 2.4. 

X J 


■.t 


t 


5830.574 mins. Since Q > d, go to 


¥e obtain, Q 
Step 3 . 

Step 3i Use branch and bound procedure. 

Set n = 1, N1 = {n} = {l}. 

For node 1, we get, R = U, 3=0, D = 0, s = 3000.00 mins. 


0 and = 1 ■¥ i. 



Calculate value of jj, from the expression given in 


Step 10 oi the algorithm. Set P. “ 



Fig. 2.2s Block-diagram for determination of 

* j * 

\x and V . . . 
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Table 2.4; Maximum Production Rate Machining 



Speed 

and Unit 

Production Time. 

__^___Index 

Group 

of 

Part 

Machining 

Speed 

ra/min. 

Unit Production Time 
min/ pc 


1 

223.6301 

6.715 

1 


2 

194.4333 

5.585 


1 

220.3509 

11.514 

2 

2 

151.0688 

12.667 


3 

151.4298 

7.164 


1 

144.3376 

8.219 

3 


2 

207.7381 

9.498 


1 

206.4631 

6.567 

4 

2 

188.0497 

11.912 


3 

112.1955 

8.338 


Step 4t 


Step 3 s 


S tep 7s 


Determine the values of the expression )/l^^ + p 




10 


¥ e'R and identify the oob for which the value 

of the expression"is minimum. For the present problem, 
J-*-» is J-ipj ond the value of expression is 6.195 min/pc. 
The maximum value of the upper bound, Z corresponds 
po node 1. Since Z 7 ^ Z, go to Step 5. 


For node 1, we find w^^ “ 50 and D = 307.25 mins. 
Since w ^2 > 

Branch the node 1 into nodes 2 and 3? and 
return to Step 4. 


(t 
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At node 2s n = n+1 = 2, N1 = {l} + {2} - {l} = {2}. 

41» '^42* ^^43^’ 


Set R *^21^ '^22.^ ^^23^ *^31^ '^32^ 


3 = 0} s = 3000.00 mins, D = 0 and Z = 0. Since R 0, 


determine {(S.+s. .)/l. . + p. for J. . e R, Identify 


ID ' ID 


ID 


iD 


the job for which the value of the expression is 

1 J 


J-i-rr is J, „ and 
ID 


minimum. For the present problem , 
value of expression is 7. 342mir/pc. Therefore, Z = 408. 

At node 3i n = n+1 = 3, N1 = {2} + {3} = {2,3}. 

Set R = {Oppj ^^21^ *^22’ ”^23’ '^31’ "^32^ '^41*' ^^42^ "^43^ * 
S = ® “ 2692.75 min, D = 307.25 min., Z = 50 


and = 0. Since R 0, determine {(S. + s...)/l.. p^ ^ 


(c) 


ID’ 


■ID 


ID 


for J. . £ R. Indentify the job for which the 

J J 


expression carries minimum value. Find that is 


for which the value of the expression is 
7.032 min/pc. Set, Z = 432. 


This time from the Step 4, we find that max Z occurs 

nsNl 

at node 2. Further branching is done from node 2 following the 


steps 4 to 7 of the algorithm. Table 2.5 gives the optimal 
results obtained. 


The optimal solution shown in Table 2.5 is compared with 
the results of the same example obtained by Hitomi and Ham [8]. 

It is observed that in their case, the total production amount 
was 356 pcs while in our case it is 372 pcs. 
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Table 2.5 j Optimal Solution for the Numerical 
Example 2. 


Index 

Group 

of Accepted 

Part 

_ Optimal 
Machining 
Speed 
m/ min . 

Unit Pro- 
duction 
Time 
min/ pc 

Quantity to 

be produced 
pcs 


1 

173.9853 

6.81 

60 

1 


2 

134.3816 

5.53 

50 

2 

3 

126.7238 

7.14 

40 


1 

117.0901 

8.06 

30 

3 


2 

134.0509 

10.08 

72 


1 

181.4794 

6.42 

40 

4 


3 

99.0188 

7.95 

80 


Total number of units produced = 372 

Total cost of production = Rs, 20,998,87 
Total time utilized = 2,999,22 mins. 


2.4.2 Case 2; Profit Maximization ii 

The problem of total profit maximization is considered 
for the following situations J 

1. The unit production time is prespecified. 

2. The machining speed which ultimately determines 
the unit production time is a decision variable. 
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3.4, 2.1 Model 3s Unit Production time is prespecifiedi 

The expression for job production time and the cons- 
traints for this problem turn out to be the same as for Model 1, 
The job production time is given by (2,12) while the various 
constraints are represented by (2.1) and (2.5). 

O bjective function s 


The function representing the total profit generated 
by the single stage production system under consideration can 
be written ast 


N. 

g 1 
Z Z 
i=l j=l 


B. . w. . X. . 
13 ij 13 


(2.42) 


where B. . is the profit per unit of job J. x. . is 0-1 
variable and is given by (2.3). The objective is to maximize 
the total profit, i.e.,Z, 


Solution 



o 

« 


In this model also, we initially find Q, i.e., the 
time required to produce all the units of all the jobs. Q is 
obtained from (2.7). IfQ<,d, all the jobs can be produced 
to their lot sizes and the optimal solution will be; 


w. 


* 


13 

7 * 


1.J V J. . e U 


= z 


max 


where , 


max 


N N. 
g 1 

Z Z B.. 1.. 
i=l j=l 


(2.43) 
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However, if Q > d, the optimal Z < Z , Under this 

^ max 

situation, we restore to B & B technique to find the optimal 
solution. 


At any stage of the B & B procedure, we select a node 

r\j 

n from all the descendents for which the upper bound value 

of the objective function, Z, is maximum. Here Z is calculated 

as the sum of the lower bound value of the objective function Z, 

and max B. . [s/{(S. l-X. + s. .)/ 1 . . + p^-^}]. Note that at 
J. .eR ^ ^ 

the fii^st node, the lower bound value Z is zero. However, at 
the subsequent nodes when a job is selected for produc- 

tion, Z is incremented by the quantity Brv wj'v for the node 
where the selected job is assigned. The B & B procedure 

branches the node n into two nodes, one in which job J'xv is 

X J 

in the set of selected jobs while in the other node Jrv is 
not in the list of selected jobs. For each of the descendents, 
we calculate the lower bound and the upper bound values of the 
objective function. The B & B procedure is terminated as soon 
as a descendant with the maximum upper bound value is obtained 
for which the upper bound and the lower bound values of the 
objective function are equal. The list of jobs in this node 
gives the optimal product mix. 


Optimizing Algorithm ; 

Following is the step by step algorithm to find the 


optimal solution. 
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Step I g 

Ste£_2^i 


Step 3 s 

Step 4 ; 

Step 5 ; 

Step 6 ; 


Calculate Q from (2.7). If Q ^ d, go to Step 2; 
otherwise, go to Step 3. 

The optimal solution isJ 
^ *• 


U 


w* . = 1 . . 

11 11 


* 


= Z 


max 


¥ J.. e S* 

11 

N 

= 2^ Z B. . 1. . 

11 11 


i=l 0=1 

Terminate the procedure. Follow branch and bound 
procedure . 

Select an initial nodes n = 1, N1 = {n} . 


For this node, set R = U, S = 0, D = 0, 


d, Z 


0 


and = 1 ¥ i. Find a 0°^ from the Qo'bs & R 

for which ~ Bw/{(Sj Xj + "*■ 

is maximum. Set Z = Btt [d/{(S:^ X-s- + sr-T;)/l-^-r + p--r} ] 
and go to Step 4. 

Find max Z and denote the corresponding node by n 
and the job in n by Jrc'. If for the node n, Z = Z, 
go to step 8; otherwise, go to Step 5. 

Find D = D + + s^ + wr-v p':^, where w^-r = max w 


11 


11 11 ' 


11 


11 ' 


w. . is varied appropriately from zero to Ipp- till 
Xo 13 

D <_ d. If w^j = 0, go to Step 6; otherwise, go to 


Step 7. 


At node n, set R = R-lJr^r}. If R = 0, set Z = Zj 

1 J •*" 

otherwise, find a oob from jobs J. . e R for which 

X J X J 
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2 = Z ^ Bj3 [s/{(SjXj sj3)/l 


IJ ' 10 


= max Y. .. Set 

J. .eR 
, 10 

+ PTilJ Sind return to Step 4. 

Step 7j Branch the node n into two nodes. 

At the first nodes Set n = n+lj N1 = N1 + {n} - {n}. 

Further, for the new node n, set R = R-{Jr~},S = S, 

10 

s = s, D = D and Z = Z, If R = 0, set Z = Z; otherwise 

find a nob J-r-? for which = max Y. .. 

J. .£R 

_ 10 
Set Z = Z + Bt-t [s/{(S^X-r + + p^-^}]. 

At the second nodes Set n = n+l, N1 = N1 + {n}. For the 


newly created node n, setR = R, S = S + s 

1 J 

D = D, Xj’ = 0 and Z = Z + wj'^ ^ij* If R = 0? set Z 
otherwise, find a nob J-r'^ for which Y-s-rr = max Y. 

^10 10 j_ _g.p 10 

Set Z = Z . BjgLs/lCSjXj . ej 3 )/lj 3 . pjj}]!' 

Return to Step 4, 

St ep 8 ; The optimal solution is given by node n. Set S''= S, 
Z = Z and stop. 

Numerical Example 3i 


rsj 

d-D, 


Consider a single stage production system in which ten 
jobs are available for productions. The jot’s are divided -into 
four groups. The basic data for the example are given in Table 2.6. 
Solution ; The solution is given using the various steps of the 
algorithm; 


Q = 3105 mins. Since Q > d, go to Step 3 



Table 2.6? Basic Data Considered for Numerical Example 
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Step 3 t Follow the branch and bound procedure. 

Set n = 1, m = {n} = {l}. 

For the newly cheated node 1, set 

R = ^12' '^21* '^22’ ^23’ '^31’ '^32* '^41* '^42’'^43^’ 

S=05 2=0, s=d, D=0 and = 0 ¥ i. Calculate 

Y. . for all J. . e R. Identify the job for v;hich Y. . 

XJ Ij IJ 

is maximum. For the present problem Y^2 maximum 
and equal to 0.9160 Rs./min. Therefore, is J^2* 

Set Z = B^2 Cd/{(S^ + ®32^/^32 ^32^^ " Rs. 2196. 00 

and go to Step 4. 

^e2_4_; max Z occurs at node 1. Since, for node 1, Z Z, 
neNl 

go to Step 5. 

Step 3 » w^2 “ pcs and p = 393 mins. Since w^2 ^ 

go to Step 7. 

Step 7 s Branch the node 1 into two nodes, 2 and 3 and return 
to Step 4. Various parameters at the nodes, 2 and 3? 
are given below s 

At Node 2i Set n = n+1 =2, N1 = {l} - {l} + {2} = {2}. 

For node n, set R = {J^l" *^12' '^21^ '^22* '^23' '^31’' '^41* 

1^2? '^43}’ S = S, Z = Oj s = 2400 mins, and D = 0. 

Since R 0, find Y. . for all J. - e R. Note that 

1 J J- J 

max Y, . is 'equal to 0.9109 Rs/min and it corres- ' 

J^jER 

poxxis to ^22* ^at 2 = ^4* ^22 ^s/l02X2*{-a22^/^22 ^22^ 

= Rs. 2178.00 



At node 3^ Set n = n+1 = 3, N1 {2} + {3} = {2,3}. 

For the newly created node n, set 

^ "^ 21 ' ‘^ 22 ’ '^ 23 ’ '^ 31 ’ ‘^ 41 " *^ 42 * '^ 43 ^’ 

S = = d ~ D = 2007 mins, D = 393 mins. 

and Z = Rs. 360.00. Since R 0, find Y. . for all 

1 J 

e R. Identify the Job for which Y. . is maximum. 

-‘-u 1 j 

max Y. . corresponds to Job J,-, . Therefore, 

J^.SR 31 

2 = Z + [s/{(S^ = Hs. 2316.0 

Now we return to Step 4 and find that max Z occurs 

neNl 

at node 3. At this node, we observe that the upper 
and the lower bound values of the objective function 
are not equal. Therefore, we go to Step 5. Steps 4 to 
of the algorithm are repeated till the optimality 
criteria is satisfied. Table 2.7 gives the optimal 
solution for this example. 

Table 2.7s Optimal Solution for the Numerical 
Example 3. 


Index_of Accepted 
Group Part 

1 


Quantity to be 
Produced 

20 


Profit 

(Rs) 


340 


2 2 50 


450 


3 30 


210 


1 

20 

240 

2 

60 

360 

1 

1 

16 

3 

40 

600 


Total profit 
Total time utilized 


Rs. 2,216.00 
2,394 mins. 
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2. 4, 2. 2 Model 4 j Machining speeds are decision variables^ 

For this model, the mathematical expressions for job 
production time, unit production time and production cost, 
and the constraints turn out to be the same as for Model 2. 

Equation (2. 14) represents the job production time while 
the constraints, are given b^r (2.1) and (2.5). The machining speeds 
corresponding to minimum-production-time and the minimum- 
production-cost are given by (2.13) and (2.16), 

Obje c tive Function . 

Our primary objective is to determine product mix and 
machining speeds such that the total profit is maximized. 
However, once the values of optimal decision variables which 
yield maximum profit are obtained, one can further manipulate 
the machining speeds such that the cost of production is minimiz 
and the maximum profit is still realized. 

The primary objective function representing the total 
profit turns out to be the same as given by (2.42). The second- 
ary criteria of minimizing the total cost of production involves 
the consideration of total production cost expression given 
hy (2.17), 

Analysis; 

In this model, the analysis is almost same as done for 
Model 2. The relationships given in (2.18) through (2.41) are 
also applicable for this model. 
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Initially, q"*", is found from (2.9) and then the time 
needed to produce all the jobs to their lot sizes at their 
maximum-production-rate machining speeds »is compared with d , 
the time available for production. Depending upon the values 
of Q and d, one of the following situations may result . 

i) £^_=,dJ 

"t 

If Q = d, then the optimal solution will be as given 

by (2.20) and the value of primary objective function Z* will 

be equal to Z 

^ max 

ii) Is 

If q"^ > d, then set p. . = p. for all jobs J. . e U, 

The optimal product mix and machine loading is obtained using 
the solution procedure of Model 3. If the slack time 6 > 0, 
then this available slack time can be advantageously utilized to 
minimize the total cost of production. In view of this, the 
modified optimal machining speeds are determined using the 
solution procedure developed for Model 2. 

iii) < d s 

" 4 - 

If Q < d, then the optimal machining speeds are 

determined corresponding to the minimum-production-cost. Thus, 

V.* = v..^®^ £ E... The values of vf. are determined 

10 10 10 11 

utilizing the slack time using the solution procedure des- 


cribed for Model 2. 





Following is the stepwise description of the algorithm: 

Find V. ^ and p. for all the jobs J. . e U from 
(2,13) and (2.12), respectively. Calculate from 
(2. 18) and (2.19). If condition, £ d, is satisfied, 
go to Step 2 I otherwise, go to Step 3. 

Step_2« The optimal solution is: 


Step 3: 



S* = U 


h - 

1 

¥ i 

•if 

no = 

1 

V J. . e S* 

1 3 

w = 

1 . 

1 

i ¥ J, . s S* 

JP ^3 

Z* = 

Ng Ni 

2 = 2 2 B. . 1. . 

®ax ^10 ^ij 


If = d, set v,*. = V. ¥ J.. e S* and stop; 

-Lj i.J 13 

otherwise, go to Step 9. 

Follow branch and bound procedure. 

Form the initial node: set n = 1, N1 = {n}. 

For the initial node n = 1, set R = U, S = 0, D = 0, 

s = d, Z = 0 and X . = 1 ¥ i. Find a job J-?— r from the 

1 13 

jobs J. . e R, for which = bW{ (S^X:5- + S7^)/l-^:3f+p 

-LJ IJ I 1 IjJ Ij 

is maximum. Set Z = B-r-^ [d/{(3-r ■t‘S^-^)/l-T-r + 

iij 1 13 ij ij 

and go to Step 4, 


11 


Find a node n from the nodes n s N1 for which Z is 

maximum. Denote the job JTT in n by Jr~. If for node n, 

1 j 1 j 

Z = Z , go to steps, else go to Step 5. 
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S^ tjep _ 6 ; 


_S^t; e_p_ 7 ; 


Step 8 ; 


Find D = D + Sp* Xp’ + spp + wpp -Q/Vrv/ (t) 
1 i ij ij 


where 


WPP = mnx vj. -5 w. . is varied from zero to Ipp till 
-1-J Ij IJ 

rv 

D ^ d. If WPP = 0, go to Step 6; otherwise, go to 

-L J 

Step 7. 

In node n, set R = R - { Jpp) o If R = 0, set Z = Zj 

10 

otherwise, find a job J-?'.’ for Vi^hich Y-s-? = max Y. 

'^10 10 . JijfR 

Set Z = Z 1 - [s/t^SjXj + 33 : 3 ) /ll 3 ■" ^in ^ 


-0 


Return to Step 4, 

Branch the node n into two nodes as given below,* 

At first nodes Set n = n+1, Ml = N1 + {n} - {n}. 

For the newly created node n, set R = R-(JPP}, S = S, 

i J 

s = s, D = D and Z = Z. If R = 0, set Z = Z; otherwise. 


find a oob J-r-^ for which Y-r-r - max Y. .. 

Set Z = Z + Bj^ [s/{(SjXj + + Pi 3 '''^^}]- 

At second nodes Set n = n+1, N1 =-Nl {n} ♦ 

For this new node n, set R = R, S = S + {Jpp}, 

-L J 


X£' = 0, s = d-D , D 


D and Z = Z + wpp Bpp, If R = 0, 


set Z = Z; otherwise, find a 0°^ which 


10 


h 

+ p 


max 


J. .eR 

(tt 


Y,,. Set Z = Z + Bt- [s/iCSjX^ + s^^)/!^- 


10 


1 1 


10 ' 


10 


10 


}]. Return to Step 4. 
The optimal solution is, 

■ S* = S 


X* 

^iO 


{ 


0 


^ J. . e S 
10 

Otherwise 
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* 


= Z 


0 


X- . { 


Ni 

If S X* . = 0 ¥ i 

j=l 


1 Otherwise 
Go to Step 9. 

Find and using (2.16) and (2.12) 

respectively. Calculate Q, 


N 


N. 

1 


Q = 2" { S X" + E (s. . + w^ . p. X. . } 

i=l a a 11 11 11 ^ 11 


If Q <_ d, set vf. = V. for J. . e S* and stop; 

■I J 1 J X o 

otherwise, go to Step 10. 

Stap^ 10 ; Select a job J. . e S*to determine initial value of p. 

1 J 


Mathematically , 

^ -]) Y . ( 


V. + V. 1/n. . 

) iO 



a 


J Using the value of u obtained from Step 10, find 

V. and p. . for all J . . e S* from (2.34) and 

11 -^11 11 

(2,12), respectively. Recalculate Q from the follow- 


ing mathematical expression; 
N 

Q = 


N. 

1 


: 2 { S . x ’. + Z ( s . . + w* . p . . ) X* . } . 

i=l 1 1 3^1 11 11 ^11 11 

If Q d, go to step 12| otherwise, stop. The present 
set of machining speeds is optimal. 
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Step 12 : Modify the value of p appropriately depending on 

whether Q<dorQ>d. IfQ<d, reduce the value 
of p and if Q > d, decrease the value of p. 

Go to Step 11. 


Numerical 


Consider the numerical example 2 alongwith the 
additional information given in Table 2.8. The objective is to 
determine optimal product mix and machining speeds to maximize 
the total profit. 


Table 2.8? Basic Datas for Numerical Example 4. 


Index 

of 

Profit per unit 

Group 

Part 

Rs./pc 


1 

10.00 

1 

2 

13.00 


1 

8.00 

2 

2 

2.00 


3 

18.00 


1 

6.50 

3 

2 

8.50 


1 

11,00 

4 

2 

14.00 


3 


5.50 
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S olution s 

Determine for all jobs J. . e U, the values of v. 

(t) 

and . These values are given in Table 2.4, 

Nov; v;ith present' set of p. = 5830.574 mins. 

t ^ 

Since Q > d, go to Step 3. 

Step , 2 ; Follow branch and bound procedure. 

Set n = 1 , N1 = {n} = {l}. 

For node 1, set R = U, S = 0 , = 1 ¥ i, Z = 0, 

s = 3000 mins, and D = 0. 


Determine for all jobs e R, the value of 

given by expression, + Pij • 

Identify the maximum value of Y. For the present 

¥• J 


problem, max Y. . corresponds to job and this is 
J. .eR 

equal to 2.225 Rs./min. Therefore, Z = B2^[d/{ (S 2 +S 22 )/^ 2 ^ 
+ ~ 6678.00. Now go to Step 4. 


Step 4s max Z is found in node 1. Its value is Rs, 6678,00. 
■ neNl 

Since Z 7 ^ Z, go to Step 5. 


Step 5 c D - 323.56 mins, and ~ ^6 pcs. Since > 0, 
go to Step 7. 


S tep 7 s Node 1 is branched into two nodes, viz, n'Ode 2 and 3. 

At node 2s Set n = n+1 = 2, N1 = {l} + { 2 } - {l}= {l}. 

Set R = R- {J 23 } = i^iiJ ^ 12 ' '^21’ '^ 22 ^ '^ 31 * '^ 32 ^ 

'^42' '^ 43 ^* S = 0, s = 3000 mins, D = 0 and Z = 0, 



Since R 0, find Y . f or ' all J, . e R. Find max Y. .. 

J. -eR 

ij 

It is observed that the maximum value of Y.^ is equal 
to 1,358 Rs./niin. corresponding to a job J]^2* 

Z - Z + [s/{(S^X^ + s^2)/1i2 ~ 

At node 3.° Set n = n+1 = 3 , N1 = {2} + {3} = {2,3}. 

For node 3, set R = R = {jqq» Jq2’ “^Rl’ '^22’ '^31’ '^32' 
"^Al^ '^42* '^43^* S = S, s = d-D = 2676,94 mins. D = D 
= 323.56 mins, and Z = 40 x 18 = Rs. 720.00, X 2 = 0, 

Since R ^ 0, calculate Y. . for all jobs J- . e R. 

i J i J 

max Y. . corresponds to job J-, ^ for which 

Yi 2 = ®12^/H2 Rs./min. 

Therefore, Z = Z + B ^2 {s/{(S 3 _ X^ + ^12^"^^^^ 

= Rs. 6375 . 00 . Now return to Step 4. 


Step 4 finds node 3 with maximum upper bound value. 
Further computations are done according to the steps 4 through l2 
of the algorithm till we get the optimal machining speeds. The 
optimal solution is listed in Table 2,9. 


Table 2.9s Optimal Solution of Numerical Example 4_. 


Index of 

Group Part 

Quantity 
to be 
produced 

Profit 

Rs. 

Optimal 

machining 

speed 

m/min. 

Optimal "unit 
production 
time 
min/ pc 


1 


' 600.00 

175.3172 

“^.787 

1 

2 

50 

650.00 

155.3731 

5.521 

2 

3 

40 

720.00 

127.4141 

7.121 


1 

19 

123.50 

117.7705 

8.034 

3 

2 

90 

765.00 

135.5620 

10.034 


1 

40 

444.00 

182.2416 

6.416 

4 

2 

50 

700.00 

159.3903 ■ 

11.225 
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Total time utilized = 2,999.66 min. 

Total profit = Rs„ 3,998.50 

Total Cost of Production = Rs. 21, 100. 00 

2.5 HEURISTIC SOLUTION PROCEDURES; 

In Sec. 2.4 we presented algorithms for single stage 
product mix and machine loading problems considering either 
the maximization of total volume of production or maximiza- 
tion of total profit. For each of these objectives two 
variations, namely, unit production time for each job is 
prespecified and the machining speeds are decision variables, 
were considered. A solution methodology for each one of these 
situations was developed. The solution methodology utilized 
the branch and bound philosoph 3 r for generating the optimal 
solution when the total production time required to produce 
all the jobs to their desired lot sizes turned out to be grea- 
ter than the available production time on the machine. In 
order to reduce the computational effort involved in the use 
of branch and bound procedure, it was felt that some heuristic 
procedure needs to be developed for identifying the jobs which 
have every high potential for inclusion in the optimal product 
mix obtained by using the algorithms given in Sec. 2.4. The 
following heuristic procedure was developed to identify the 
set of jobs having high potential for inclusion in the final 
solution. 



As soon as ishe first step of any of the algorithms 
given in Section 2.4 reveals that the total available time 
on the machine is less than the total time required to produce 
all the jobs to their lot sizes, contribution of each job 
towards primary objective function, e.g,, maximization of total 
amount of production or the total profit is calculated. The 
job with the highest contribution is identified and selected 
for inclusion is the list of selected parts provided the time 
required to produce the job to its lot size is less than or 
equal to the available production-time. After the first lot 
is selected, the lot with the next highest contribution is 
identified and the procedure is repeated till no job can be 
selected for inclusion in this list of selected jobs. 

If the number of selected jobs in this list turns out 
to be less than or equal to two, this list is made empty and 
the branch and bound procedure, i.e, step 3 onwards of the 
various algorithms are followed as such. However, if the list 
of selected jobs contains more than two jobs, the last two 
jobs are removed from the list. The Step 3 of the branch and 
bound procedure is suitably modified for each of the algori- 
thms given in Section 2.4 to incorporate the heuristic 
procedure into the branch and bound approach. The remaining 
steps of the branch and bound procedures, i.e., Step 4 
onwards remain the same as in the original algorithms. The 
modified Step 3 for each case is presented below. 
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2 « 5 . 1 .1 *„ of Tolsal Unit ; s Produced : 

2. 5 • 1*1 Coiso l(a); Prespecified Unit Production Time 

Hp;uri_stic^ A?-£. 9 rithra j 

The Stop 3 of the original algorithm given in Sec. 
2. A, 1,1 in modified as follows j 


SJep, _3 ..1 . 


5.2: 




Step 


Select an initial node? n = 1, N1 = {n} . For this 

nodes set R=U, S=0, 0=0, s=d 5 Z=0 and 

= 1 V i . Go to Step 3.2. 

Find from nobs J. . c R, a nob J-.-t for which 

[(St Xt + st - t )/ 1 t -» + Ptt} is minimum. Find 
^ 1 1 10 10 10 

'D = D + St Xt + Stt h- Wtt Ptt , v/herG Wtt w .jW . 

11 10 10 10 10 10 10 

varied from zero to Irr till D <_ d. Go to Step 3*3. 

i J 




set S = S + {Jtt} 

1. . and go to Step 8. If wt? = 0 or 


If Wtt = It'S end D = dj ocu i_i — u t 

-i- V ***“ tj 


iO 


S + {JFj}. Xt =0, D = DsS 


d»D 


10 

vj^T-r < i-T , go to step 3.^5 otherwises set 
10 iO 

R = R _ {j_}, s 

and go to Step 3.2. 

If S = 0 and Wtt = 0, then set Z = 0, S*= 0 and stop. 

JL J 

If set S contains more than two JoPs, remove the 

last two oofes selected from the list S and include 

these jobs in R, otherwise, make S empty. Set 

X. = 1 for those groups G. ’of which no job is 
i 

included in S. Set Z = 2 1. D = { 2 

J. xS J.. xS 


N, 


+ w . p. •) -i- 2® S.X.} and s = d~D, Find a job 

iO 10 i=i 1 1 


10 


10 


is 
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Jj- from jobs J.^ e R for which {(S^X^ + 
is minimum. Set Z = Z + [s/{(S^X-^ + s--r)/l-^ + p^-^}3 
and go to Step 4. 

2.5.1«2 Case l(b)i Machining Speeds are Decision Variables^ 
He uris tic A lgorithm ; 

Step 3 remains the same as given in Sec. 2. 5 ‘I.! 

except the fact that instead of p. given by (2.18) 

J 1 J 

is used. 

2.5.2 Ca se 2; Maximization of Total Profit ; 

2. 5. 2.1 Case 2(a) i' Prespecified Unit Production Times 
HeurisJ: ic Al gori thm s 

The Step 3 of the optimizing algorithm given in 
Sec. 2. 4. 2.1 is modified as followss 
St ep 3 .1s 

Step ,3r2.° 


Step 3.3 s 


Select initial list, set n = 1, N1 = {n}. 

For node 1, set R = U, S = 0, D = 0, s = d, 2 = 0 


and X^ = 1 ¥ i. Go to . Step 3.2. 


Find a job Jrr from J. . e R for which B*^/{(S^XT+S:^-v)/ljt^ 

1 j 


10 


tj 


1 i 10' 


'10 


D + S-X^ + S-- 
11 lO 

is varied from 


+ p-r-r} is maximum. Calculate D 
10 ^ 

+ pjj, where Wj^ = max w^^, w^^ 

zero to Irr till H < d. Go to Step 3.3. 

J 

If ■w-r-r = It-^ and D = d, set S = S + 

lO 10 ■^0 

Z - s B. . 1. . and go to Step 8. If = 0 or 

- J eS 

•W_, < 1-,, go to step 3.4; otherwise, set R = 

it] i 0 
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3 o '(• Xj - 0, D = D, s = d-D and go to 

Step 3.2. 


Step_3^,. If S = 0 and = 0, set 2^^= 0, s''= 0 and stop. If 

the set S contains more than two jot), remove the last 
tv;o jobs put in list S and include them in Rj 
otherwise, set the list 3=0. = 1 for those 

group of which no job is included in S. 


Set Z = 2 w. .B..,D={2 

J..eS J. .eS 

10 10 


(s. . + w. . p. .) 
13 10 ^10 


N. 


i=l 


. X. } and s = d-D. Find a job J-r-r from the jobs 
L 1 ^10 


J. 


iO 


e R for which Bjy{(SjXj ^ Pjg} is 


maximum. 


Set 2 = 2+ [s/{(S^X^ + Sj^)/l-, - p,-.}] 


1 1 


10 10 


10 ' 


and go to Step 4. 


2. 5. 2.1 Case 2(b ) j Machining speeds are decision variables. 
Heur istic Algorithm s 


Step 3 remains the same as that given in Sec. 2.5.2. 1 
except the fact that instead of p. •» p. given by (2.16) 

1 J 1 J 

is used. 


2.5.3 C ompu t at i ona 1 Exp e r i enc e s 

To check suitability and computational efficiency of 
the heuristic procedure, numerical examples given in Sec. 2.4. 
were solved using heuristic algorithms. It was observed that 
in all the cases, the solutions obtained by heuristic proce- 
dures and the optimizing solution procedures were the same. 
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Table 2.10 gives the CPU time taken on DEC-1090 computer system 
for the various examples when solved using the optimizing 
solution procedures and the heuristic approaches. It is 
observed that for all the examples considered the CPU time 
requirement was lower for the heuristic approach as compared 
to the optimizing procedure used for solving the example. 

Table 2.10; CPU Tlm.e Taken by DEC-1090 Computer 
System. 

Problem 

Optimizing Solution Heuristic 

Procedure Procedure 

Numerical 

Example 1 0.13 0.10 

Numerical 

Example 2 0.16 0.14 

Numerical 

Example 3 0.13 0.11 

Numerical 

Example 4 \ 0.20 0.16 

2.5.4 An Illustration Demonstrating Use of Heuristic 
i^o'cedure i 

Numerical example 1 presented in Sec. 2.4.1. 1 is 
considered to demonstrate the heuristic solution procedure. 
Solu tion s 

The following is the step-wise solution of the example 
using the heuristic approach. 
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^31^ '^32’ 


St ep, Is Q = 3105 mins, and d = 2400 mins (given). 

Since Q > d, go to Step 3. 

stop 3 . I s Select an initial nodes set n = 1, N1 = {n} = {l}. 

For node 1, set R = (JxL, '^21' '^22’ ^2' 

J^l, '^ 43 ^’ S = 0, D = 0, s = 2400 mins, 2 = 0 

and = 1 ¥ i, go to Step 3.2. 

Step 3. 2!. Calculate min {(S.X. -i- s. .)/!.• + P,--;}* 2ts value 
" J^^eR ^ ^ 10 11 10 

is 6.55 min/pc corresponding to job 1 ^ 2 * 

that D = 393.00 mins for ■w ^2 “ ^6 pcs. Go to Step 3.3' 

Sto p 3 . 3s Since D < d and = ^ 32 ' ^ "" ^-{^ 32 ^ ’ 




0 J > l/i O > * 


' 12 ’ '" 21 ’ ^" 23 ’ "" 31 ’ "" 41 ’ “^ 42 ’ ^ 43 - 
S = S + (1^2^ = ^"^32^’ Xj = 0 and D = D = 393.00 mins. 

Go to Step 3 . 2 . 

Now we recalculate {(Sj^X^^ + Sj^j)/lj^j + ^ij^ 

J. . e R. It is observed that the minimum value of the 
3. ^ 

expression is 6.8 min/pc corresponding to o®!* '^12* 

We find that 'W’ 2^2 ~ pcs and D = 461.00 mins. 

Go to Step 3 . 3 . 


Since '''^]_2 ~ ^12 D < d, set R 


{Jll, J2;j_, J22»J23» 


12 ’ 


1 ^ 2 } » = 6 and 


■^ 31 ’ ‘^ 41 ’ ^^ 42 ’ *^ 43 ^’ ^ 

D = 461.00 mins. 

Similarly following Step 3.2 and 3.3 iteratively, we 
find that, R = {Joi* '^41’ ^^42^’ ^ ^'^11’ '^12’ '^22’ '^23 


J 


^ xiv^O^^i, D = 2289.00 mins, and 

31 ’ '^ 32 ’ ^43’’ ■' 


w 


’42 


= 10. Since w ^2 ^42’ 
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'ree-Diagram obtained for Numerical Example 1 using Heuristic 
■pproach. 
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S tep 3«^ « We remove the last two jobs and from S and 

include them in R. Now set R = {J22’ '^31' “^hl' "^42 *'^43^ 

S = '^22* ‘^31’ '^32^' Z = 180^ 

D = 1379.00 mins, and s = 1021.00 mins. 


Find min {(S.X. + s. .)/l. . + Zts value is 

J sR ^ ^ 

IL Q 

12.25 min/pc corresponding to job Set 

Z = 2 + [ s/iCS^X^ + + P3 i}3 = 263 and 

go to Step 4. 


Now we follow branch and boimd procedure iterating 
steps 4 through 7 as discussed in Sec.. 2. 4. 1.1. The final 
product mix obtained is as shown in Table 2.2. Fertile purpose of 
illustration a portion the tree-diagram obtained by the use of 
heuristic algorithm is depicted in Fig ^-2.3. 



CHAPTER III 


MULTISTAGE PROBLEM 

3.1 STATEMENT OF THE PROBLEMS 

Consider a multi-stage production system v;ith Ns 

production stages. Each production stage has limited 

available processing time, N ,t 3 rpes of jobs are available 

for production. Based on the part-family concept of group 

technology, the jobs have been classified into N groups. 

S 

The i-th group,, Gj^ (i = 1,2,..., N ) has N. jobs. Each job 

J. . (j = 1,2,..., N. j i = 1,2,..., N_) has o lot size 

1 J 1 s 

associated with it. 

The objective is to optimally select the jobs and 
determine the quantities of each to be produced such that 
the desired criteria of optimization is satisfied subject to 
the restriction that the total processing time available on 
each stage is not exceeded. The optimization criteria is 
either to maximize the number of units produced or maximize 
the total profit. 

The unit processing time required to process a 
single unit of the job at a given stage is either pre speci- 
fied or determined based on optimal selection of machining 
speed., The optimal selection of machining speed is based on 
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minimum cost of production, 

3.2 ASSUMPTIONS J 

The various assumptions are; 

i) The jobs to be processed on a multi-stage production 
system d" not have any prespecified sequence of c-pei’ations. 
Thus, job operatiions do net have any technological ordering. 

ii) The various jobs do not necessarily require 
pi’ocessing at all stages of the production system. 

iii) Jobs do not have to be necessarily pr^^duced to 
thoir lot sizes. 

iv) Once e job is undertaken for production on a stage, 
all the units of this job are completed first before taking 
up the next job. 

v) Each unit of the job must be processed to completion. 

vi) Each jjroduction stage has only one machine. 

vii) Only one unit of the job can be processed on the 
machine at a time. 

viii) A limited amount of processing time is available on 
each machine (each stage). The processing time is the actual 
productive time available on the machine excluding the time 
lost due to break-downs, etc, 

3.3 NOMENCLilTURE: 

Following are the notations used for the development 
of the mathematical models* 



^ijk ~ time for job J. . at k-th stage 

1 3 

(min/pc) 

^ijk replacement time for job at k-th 

stage (min/pc). 

B^. = profit per unit of job J. . (Rs./pc). 

*^0 

= 1-minute tool life machining speed for job 
Jfj at k-th stage (m/min) 

CUk = total utilized processing tim.e of the machine 
at k-th stage (min) 

dk = available processing time at k-th srage (min) 
Eijk = high-efficiency machining speed range for job 
J . . at k-th stage (m/min) 

J- J 

= i-th group or part-family 

i = group index (i = 1,2,..., N ) 

S 

j = job index (j = 1,2,,.., N^^) 




N 

“g 

h 

Ns 


P 


ijk 


= j-th job in i-th group 

= index of the production stage (k = 1,2,..,, Ns) 
= lot size of job J. . (pcs) 

-i- J 

= slope constant of the Taylor tool- life equation 
for k-th production stage 
= total number of jobs available for production 
s= total number of part-families or groups 
= total number of jobs in group 

= total niomber of stages in the production system 
= unit processing time for at k-th stage 
(min/pc) 
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(c) . . 

Pijk ~ roii^roum-producti on-cost unit processing 
(t) 


time for J. . at k-th stage (min/ pc) 
1 J 




maximum-production-rate unit processing 
time for J. . at k-th stage (min/pc) 

-L J 


Pij(Vi^)= machining speed dependent unit processing 


^ijk 

S., 

ik 

^iak 

'^iik 

U 


time for J. . at k-th stage (min/pc) 

= 30b set-up time at k-th stage for J. - (min) 

X o 

= group set-up time at k-th stage for group 
(min) 

= actual machining time at k-th stage for 
(min/pc) 

= tool life at k-th stage for J. . (min/ edge) 

-L ^ 

= set of all parts available for production 


^ijk 


(c) 




(e) 




(t) 


w. 


11 


X 


13 


X.. 


Y 


{Ji^ 3 - 1 , 2 ,..., N^i i - 1 , 2 ,..., Ng) 
machining speed for J . . at k-th stage (m/min) 

13 

minimum-production- cost machining speed for 
J. . at k-th stage (m/min) 

13 

high efficiency machining speed for at 

k-th stage (m/min) 

maximum-production-rate machining speed 
used at k-th stage for J. . (m/min) 

X J 

quantity of J. . to be produced (pcs) 

J- ^ 

0-1 type variables for 

0-1 type variable for 

total cost of production (Rs. ) 
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Y ’ 

Z 

z 

z 


max 


a 


k 


*^ijk 


Y 


ijk 


idk 

0 


1 ^ 


-X- 


= total variable cost of production (Rs.) 

= value of objective function 
= lower bound on objective function value 
upper bound on objective function value 
= maximum value of objective function 
= direct labour and overhead cost at k— th 
stage (Rs./min) 

= machining overhead for at k-th stage 
(Rs./min) 

= cost of tool at k-th stage used for 
(Rs./ edge) 

= slack time at k-th stage (min) 

= machining speed constant for at k-th stage 

= null set 

= Lagrange multiplier 
= shows optimal value 


Following additional notations will be used in the 
development of solution procedure: 

= available processing time utilized at k-th 
stage (min) 

n = all nodes belonging to N1 

N1 = set of all descendants 

R = set of remaining job at any node 

Sj^ = time still available for processing at 

k-th stage (min) 
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S = set of selected jobs at any node 

= shows a candidate for acceptance or rejection 
[a] = is a Gaussian notation implying a greatest 
integer less than or equal to A. 


3.4 matpiematic/il models and solution procedures j 

Mathematical models for the multistage problem stated 
in Sec. 3*1 presented along with optimal solution methodo- 
logies. The maximizaticn of profit is considered as the cri- 
teria f jr optimization. Separate models and solution methodo- 
. logics are presented for the following two cases s 

i) The unit processing time is prespecified, 

ii) The machining speed which ultimately determines 
the unit processing tim-e is a decision variable. 

3.4.1 Total Profit Maximization When Unit Processing Time_is_, 
Pre specified ? 

3. 4. 1.1 Constraints? 

The following constraints are considered? 


i ) nt.il-ization constraint ? 

This constraint ensures that the total time needed 
to process the various jobs at any stage does not exceed the 
available capacity of that production stage. Mathematically, it 

can be written as? 


g 

i=l 


(S., X. + ^ 

IK 1 0 


N. 

I 


JL p 


. X. . ) < 

Ic == 1, 2, » • 


Ns 


(3.1) 
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whore , 


processing tine for processing 


units of job J . . at the k-th stage and is given by, 

X J 




(5.2) 




jk " °ijk ■'* ^ij Pijk ’ 

1 if job J^. is selected for production 


{ 


and, 


ID 

0 if job is not produced 


N. 


if 2 X. ■ > 1 




1=1 

N. 


il 


(5.3) 


(3.4) 


0 


if E X. . = 0 

j=l 


ii) Lo t Size Constraint ? 


The quantity w^^ of job produced should net 


exceed the lot size specified for it. Thus, 


0 < Wy < lij 


J . . a d 
11 


(3.5) 


3. 4. 1,2 Objective functip_rj_ 


The total profit can be expressed ass 
Z 


Z j: B, . 

1=1 j=i ^ 


(3.6) 


where, Is the profit per unit of job is a 0-1 

variable defined by (3.3). The objective is to maxinilzo the 

total profit, i.e. ,Z. 
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3. 4. 1.3 Solution procedure s 


Lot represents the time needed to process all 
the jobs to their lot sizes at the h-th sta[';e of production. 
Then, 




‘k 


Ng 

^ '^ijk ^ij 

k = 1, 2, . . . , Ns 


(3.7) 


If for all the production stages the condition 

jr 

Q <_ d is satisfied, then all the jobs will be produced to 
their lot sizes. Thus, the optimal solution will bes 


where, 


ha 

n*a 

* 

Z 


= 1 

= 1 

= z. 


Ij 


¥ J. . e U 
13 

¥ J. . e U 
13 


max 


N. 


(3.8) 

(3.9) 

(3.10) 


max 


i=l ,1=1 


N. 

1 

^ ha 


If the condition is not satisfied for any 

of the production stages, then any one of the following 
situations may arises 

i) All the jobs except one of them, are produced to 


their lot sizes. 

ii) All the jobs are not selected for production. 

In this situation two cases are possibles 
a) Jobs selected for production are produced to 


their lot size. 
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b) One of the selected ;3obs can not be produced 
to its lot size. 

In all the above situations, it is obvious that 
^ ^ ^max* branch and bound (B & B) procedure is developed 

to determine its optimal product mix and machine loading for 
each of these cases. 


The B & B procedure is used to branch a node n for 
which the upper bound value Z, is maximum. For each of the 
descendent obtained after branching a node, the upper botind 
value Z and the lower bound value Z is determined. To obtain 
Z, a ^ob from the remaining jobs is identified for which 


B. . N. . is maximum. The corresponding job is denoted as 
Mathematically, N. . is expressed ass 



min [ 
k 





( 3 . 11 ) 


where s^^ is the time available for production at that node. 

The job J'TV is included in the list of selected jobs in one 

^ J 

of the two nodes obtained from branching of a node having 


maximum upper bound value. For the descendent, for which the 
job Jrr is included in the list of selected jobs, the lower 
bound value Z, which at the beginning of B & B procedure was 
set equal to zero, is incremented by B^ w^. The list of the 
selected parts corresponding to the node for which the upper 
bound value is maximum among all the nodes and equal to the 
lower bound value at the same node, gives the optimal product 
mix. At this point the B & B procedure is fathomed. 



3. 4, 1.4 ^timiziriF: algorithm ; 

The various steps of the algorithm ares 
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Step 1 

Step 2 

Step 3 


Step 4 

Step 3 


Find Q^, for each stage from (3.7) . If ^ ^ k, 

go to Step 2j otherwise, go to Step 3, 

Produce all the jobs to their lot sizes and terminate 
the procedure. The optimal solution is given by the 
fallowing relationships s 


X 


id 


^ * 
1, V/. . 

i:i 


1 . . ¥ J . e U and Z ' = Z 

13 13 max 


Follow branch and bound procedures 

Select .an initial node? set n = 1, N1 = {n}. 

For this node, sot R = U, S = 0, 

Find, 

S-, 


= 0 ¥ k, Z 


0 and = 1 ¥ i, 


Find N . = min [ 

fik fij ®ijk 

"■ ifr 


] ¥ J. . s R- 


+ p 


ijk 


Bind max g denote the corresponding job 

JijSR 13 13 

by Jrr. Set, ^ = Btp Nr*: and go to Step 4. 

*13 d 3 13 

Find a node from the nodes n e N1 for which Z is 
maximm. Denote this node by n and in n by 

For the node n, if Z = Z go to Step 8? otherwise, go 
to Step 5. 

Find for each stage \ “ Dk + + "Oj Pljk' 


max w. w. . is varied from zero to 1'^ 

1 J X J J 


where w^^ 

till < d^ ¥ k. If w^ = 0, go to Step 6? otherwise, 
go to Step 7. 
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Ste p 6 ; 


St ep 7 5 


Step 8 s 


Remove the job from the set R of the node n. 

Sot R = R - If R = 0» set 2 = Z| otherwisej 

find, ®ij ^ij* denote the corresponding job 

by and set Z = 2 + %j* Return to Step 4. 

Branch the node n into two nodes, as given belcws 
At first Nodes 


Set n = n+1, N1 ■= N1 + {n} — {n}. 

For the newly ere ted node n, set R = R - { , 


S = 0, - Sj^ and 

and Z = ^. If R = 0, set 2 


¥ k, N. . = h. . ¥ J. . e R 
IS 10 10 

= Z| otherwise, find 


max B. . N. denote the corresponding job by J-r-r 

J. .eR 
10 „ 

and set 2 = 2 + Bt-^ N-!-r, 

” 10 10 

At second Nodes 

Set n = n+1 and N1 = N1 + {n}. 

For the newly created node n, set R = R, S = St 

J 


“ ^k ■ °k ^ 


X. =0 and 2 = 2 t 


+ Bo^ w-TY. If R = 0, set 2 = 2? otherwise, recalcu- 
10 10 

late the value of N. . ¥ J . . e R and find max B N . . 

Xj ij IJ IJ 

Denote the corresponding job by Jtt and 

IJ 

sot Z — ^ 

X J X 

Return to Step 4. 

The optimal solution is' contained in node n. 


* 


Set, S = S, 2 


2 and stop . 


The value of wo^ required in Step 5 is determined using 
X 

the following relationship. 
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w~~ = 

min [ 
k 

itihL. 

- Sr, 

XK 

13 


^ijk 


If, wrv 
10 

\ 1 

iO' 

we set 

id 


3. 4. 1.5 Numerical Ex.- u nple 5 ; 

Consider a two stage production system for which ten 
types of jobs arc available for production. These jobs are classi- 
fied into iour groups. The basic data for the problem is given 
in Table 3.1. The objective is to determine optimal machine 
loading and product mix. 

3. 4. 1.6 Solution j 


Using the various steps of the algorithm, we obtain, 

1 ^ = 3105.00 mins, and Q2 = 2995.00 mins. 

Since > d^ ¥ k, go to Step 3. 

St op 2 ! Select an initial nodes Set n = 1, N1 = {n} = {l}. 

Foi node 1, set R = ^'^ll* '^ 12 ^ '^ 21 ^ '^ 22 ^ *^23^ '^32’ 

'^42’ '^43^ » S = 0, = 2400.00 mins, = 2400.00 mins, 

Z = 0 and X. = 1 ¥ i. Calculate K. . ¥ J. . e R, Find 

*** X Xj XQ 

max B. . N. .. Its value is Rs, 4080.00 corresponding to 
J. .cR 

job Ujo* Therefore, Z = B-,2 N^2 “ ^080,00. 

Go to Step 4. 

Stop 4 s Since for node 1, Z Z, go to Step 5. 

Step 5: For w^2 ” “ 68.00 mins and D2 = 94.00 mins. 

Since < dj^ ¥ k, go to Step 7. 

Step 7 s Branch node 1 into two nodes, viz, nodes 2 and 3. 

At node 2s Set n = n+1 = 2, N1 = {l} + {2} -- {l} = {2}. 
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Tabic 3.2s Optimal Solution of Numerical Example 5. 


Index of 

Group 

Accepted 

Job 

Quantity to 
’be produced 

Profit 

obtained 

Rs . 


1 

30 

195.00 

1 





2 

10 

160.00 


1 

20 

390.00 

2 

2 

50 

625 . 00 


3 

30 

270.00 

3 

2 

60 

840.00 

4 

3 

40 

592.00 

Total 

number of 

units produced = 240 

Total profit obtained 

= Rs. 3,072.00 

Total 

processing time utilized: 


At stage 1 

= 2,394.00 minso 


At stage 2 

= 2,386.00 mins. 
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and '6=0, We find that .s^ = 82 = 2400,00 mins, 

Z = 0 and = 0 ¥ i. Since R ^ 0, find 

max B. . N. .. Its value is Rs. 3025.00 correspond- 
^0 10 

ing to a job J22* Set Z = Rs. 3025.00. 

At node 3s n = n+1 = 3, N1 = {2} + {3} = {2,3}. 


R 


S 


J22, ^22’ ^23’ "^31’ ^^32’ '^41’ '^42’ '^43^ 
{J 12 }. We find that = 2332.00 mins, 


and 


D2 = 2306.00 min, = 68.00 mins., = 94.00 mins, 

Z = Rs. 160.00 and = 0. Recalculate the value of 
W. . ^ £ R* Find max g ^ value jcorres- 

ponds to a job J22* Therefore, ^ = Z + B22 ^22 ~ 


- Rs. 3110.0. 


Return to Step 4. 

From Step 4, we find that max Z occurs at node 3. 

nehl 

Therof'ire, the further branching takes place at this node. 

The steps 4 to 7 of the algorithm are repeated till the 
branching procedure is fathomed. The optimal solution obtained 

is given in Table 3.2. 

3.4.2 Total Profit Maximization wh en Maohinlng Speeds 
are Decision Variables . 

Before developing a mathematical model and solution 
procedure for the problem, we develop mathematical expressions 
far jo-b processing time and cost for processing a single unit 


at a given stage of production. 



3 . ^ 1 Job processing tim e; 

The unit processing timej Pj_j]^ ^^ijk^’ required to 
prucosG a single unit of the job at the k-th production 
stage' is rx function of machining speed for that job at 
k-th prolAiction stage. Mathematically 9 


•'ijk ^"^ijk^ 


+ ^i.1k 

a. + t. + b. fn-"^ 
^ijk Ilk Ilk 


Mr. 


^1/n. -1) 


I ±/ ii . -T, 

hua hi 

(5.12) 


ndk 

¥ J.j e U 


whcT,.., b.ji^ and respectively represent the 

pro;:.- '.rail'll time, tuol replacement time, and actual oaclilnln£ 

time on t-th production stage. ‘^ijk 

parameters of the Taylor tool life equation. 

l,et V. . represent the machining speed at which 

P, Cv...p Is'mLimum. is obtained by setting the 

first dirtvative of with respect to v^j, equal 

to zor *. 


(t) _ r i = 1 » 2 , 

"^ijk “ ?^ijk^ 


« 4> « 9 


(3.13) 


, % 


3 “ 


■+- -r -irih T are processed on k-th stage at a speed 
If w. . units nf 10b are pxoi. 

X 3 ■ ‘ 


V 


ijk 


the job processing time, is given by, 


\ * 

ijk 


®ljk ■'■ na ^ijk ‘^''ijh^ 


(3.1^) 
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• 3 . 4 , 2 . 2 Production Coat. ■ 


The cost of processing one unit of job J. . at stage k, 
^ijk also a function of the machining speed v. 

at th<3t stages Mathematically, it is expressed as; 

* '“k hjk Ujk.) ifeu- 

^ijk 

(3.15) 

whore Yijk’ respectively represent the direct, labour 

cost, machining overhead cost and tool cost for job J. . at 

i<D 

pr^jducti<in stage k. 

(c) 

Let ^ represent the ’Minimum- Production-Cost 

Machining Speed’. It is obtained by setting the first deri- 
vative ef ( 3 . 15 ) with respect to equal to zero. Mathe- 

matically, 


^ijk 


C.,v { 


^°-k h.ik^ j' 

^“k hjk Ujk^ 

(3.16) 


3 . 4*2. 3 Ob.iectr 


Our primary objective is to determine product mix 
and machining speeds such that the total profit is maximized. 
However, once the values of optimal decision variables which 
yield maximum profit are obtained, the machining speeds can 
be further manipulated such that fhe cost of production is 
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I 


minimized without sacrificing the maximum profit. The primary 
oboective function turns out to be the same as given by (3.6). 
The secondary criteria of minimizing the total cast of produc- 
tion involves the total production cost expression which is 
given as: 


Y = 


Ns Ng 

s 

k=l i 


N. 

1 


^“k hk T (“k ®iok 


+ w.. C }..1 VV..1 

10 ^lOk ^ lok 




(3.17) 


3 . ^ . 2 . 4 ^aly si s 


I'or all the oo’^s, the minimum-producti on-time machining 
speed and unit processing time at the various production stages 
are calculated from (3.13) and (3.12) , respectively. The follow- 
ing expression is used to calculate the jot processing time for 
each Job at each stage considering the fact that the each unit 
in the lot associated with that Jot is produced at minimum unit 
processing time. 


P = S ' + 1 O S (v (t)") 

ijk ^ijk ij ^ij^ ^^ijk ^ 


(3.18) 


X 

Mathematically, , the time needed by k-th production stage 
to process all the Jobs to their lot sizes at minimum-production- 
time machining speeds, can be expressed as, 

N 


N. 

< “ ^ hjk 


(t) 


), k = 1,2,..., Ns (3.19) 


Depending upon the value of , there are three 
possible situations s V k, > dj^ for at least 
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one of the production stages and Q, ”*" < d, k. These three 

K K 

cases are discussed separately, 

= dj^ ¥• k (k = 1,2, ♦,.j Ns), then produce all 
the jobs to their lot sizes (at the minimum-producti on-time 
machining speeds). In this situation, the optimal product 
mix and the cutting parameters are given by the following 
relationships; 

^ij " ^ 

w. * = 1. . r V J. . e U (3.20) 

10 IJ 10 

and "'^ijk ~ "^ijk ^ ~ 1,2,,.., Ns 

ii) > dy. for at least one of the prod uc tion stages ; 

In this case, the processing time available on at least 
one of the production stages is not sufficient to process all 
the jobs to their lot sizes. In fact, any one of the following 
situations may result; 

i) All the jot's, except one are produced to their lot sizes. 

ii) All the jobs are not produced and 

a) the jobs selected are produced to their lot sizes. 

b) all these jobs except one are produced to their 
lot sizes. 

In all the above situations, the optimal product mix 

be obtained using the following approach. 


and machining speeds can 
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First of all, we determine the optimal machine loading 
Lind product mix based on the minimum“production~time speeds at 
c.nch stage using the solution methodology discussed in Sec.3»^»l* 
Lat and £ S ) respectively represent the set of 

jidas accepted for production and the quantities of each Job to 
bf’ r)roduced. Following is the optimal solution. 


x.i 


{ 


R-irther, 


^i 


* 


= { 


1 J. . £ 

0 otherwise 


% 

1 If 2 X. . > 1 

0=1 

0 Otherwise 


( 3 . 21 ) 


( 3 . 22 ) 


I'Vtr c-Mch production stage, the total time consumed Cu^^ and 

th*.' slack time 6^ are calculated from the following expressions. 


Cu, 


Ni 


and 


i=l 


(t) 




Cu 


k 


(3.23) 

(3.24) 


If = 0, then all the selected Jobs are processed at their 
miniuum-production-timo machining speeds (at that stage). 
Obviously, 

(■t) J1 .T . TT (3.25) 


* 


= V. 


¥ J. . £ U 
10 


ijk “iOk 

If 6 > 0, it implies that there is some slack time 

available at the k-th stage. This slack time can be advantage- 
ously utilized to minimize the total cost of production. Since 


"'^i Jk 


(t) 


> v^ and q. > ^Jk ^^ijk'"')’ 


(t)^ 


(c) 
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« 

rcjasf^nablG to assiime that the cost of processing a single 

unit of a job is least at the minimum-producti on-cost machining 

spc'Cfl and it increases with an increase in the machining speed 

the minimum-production-time machining speed. Therefore, 

the ulack time can be utilized in lowering down the cost of 

producticni by decreasing the machining speed from 

(c) _ _ * „ ("t) 


V, 


Thus, the objective is to determine » 


"^ijk 


such that the cost of production is minimized. 


ijk 

< V.'-v 

i K. 

The problem is formulated as a nonlinear program as given belowi 

N. 


10 


N, 

Ns g * / 

Minimize Y = 1 2 ®ijk 

k=l i=l 0=1 




( 3 . 26 ) 


GUbjoCt to 
N 


g 


thk h 

1=1 


* 


It 

^®ijk ^ij 


?ijk ^^ijk^'' ^ij ^ ~ 


Ns 


(3.27) 


k — lySj ***. 

There are certain constant terms in (3*26) and (3*27) 
which can be ignored for the purpose of optimization. The 
problem can be restated as: 


Minimize Y ' 


N N. 

Nsgi^ I \ * i-x 

^"idk) njk (3 

k=l 1=1 0=1 


subject to, 
N_ 


N. 


i=l 0=1 


* 

"13 P 


4 


ijk ^^ijk^ ^iO 


< do, k = 1,2, . 


.28) 


Ns 

(3.29) 
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where, 

Y’ 


Ns N^ 

“ " k=l ifi hk h* - 4^ % 3,^, X*. ) (3.30) 


and, 


N N. 

1 f ^ f -:^ 1 

k = d,- L (S,,X. +Z s. X* 'S k - T o 

^ 1=1 ^ i.ik 1 .1 ^ K ^ 1,29 ♦o.y Ns ( 3 # 31 ) 


0=1 -10 


For the problem represented by (3.28) and (3.29), a Lagran- 
gian function of the following form is obtained. 


'^nak- '^k) 


Ns ^g Ns 

= 2 E Z w.*.'a (v '1 X * ■ 

k=l i=l 0=1 ^ iik'' "^iOk 

NfiJ \ N. 

+ 2 ) T, { 2 S V / . . p , . . ( V . . , ) X r . 

k=l i=l j=l ijk ij 

+ - do } 

v/hc.re <'...nd represent the Lagrangian multiplier and slack 
variable, respectively. For optimal values of k ^ ’ 

we set = 0. It gives 


'^hk - 


0 


iOk 


or. 


c 


+ p, 


^ ’-^Idk * K hsk *^ljk) 


^iOk 

O IJ 

^id 
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^k ^ vT^ 


i0k_4ok 
i/n. ., 
r ijk 

idk 


^ (l/n^j]^-2) 


) = 0 


From 'tho above equation, we get 


/ , 1/n. 

(1/n. -1) Y . (v* /C ) 

■/ ink ^ '^i.i ^^i.ik / ink-^ 


■1 Y b. (v*.^ /C. OK 


1,2,..., Ns (3.32) 


^ijk 


C. { 
lOk ^ 


“k '^k 




^k ^ijk ■*'^iok^^^/^i3k“^^ 


i = 1,2, . • ♦ c N 


3 ^ 1,2,*64I, 

Ic ~ 1,2, «ro*, Ns 


(3.33) 


SI milnrly 


, setting = 0 yields. 


iM 

g N. 

no Piok <nok ) k*. . k = 1,2,..., No 

(3.34) 


an;i 3 z/ 0K, = 0 gives, 


— 0 ]rc -r.^1, 2,«*«, Ns 


:3.35) 


If we substitute v* , = v. in (3.32), we find that 

= 0. Similarly, if we substitute = "^ijk (3.32), 

we obtain = <». Hence, it is obvious that, 


^ ^ ° uok " <nok'"^- nok^'" 3 . k ■= 1.2 ,... ns 


(3.36) 
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Since is non-negative, it guarantees the minimization of 

cost of processing selected jobs on k-th production stage. 

If 4 = 0, then from (3.33), = v. ^ j ^ g* . 

-L J 

Hc)Wc-vor, if > 0, then the optimal values of and 

e ^ijk each of the production stage is deteimiined 

frnm ( 3 . 32 ) and (3,33) , respectively, by using hit and trial 
aipproach. E. is the high-efficiency machining speed range for 
j<.»b J.. ■ at stage k. Mathematically, 

IL 3 


E. = [v. 
^i 3 k 13 k 


(t) -1 
V. M J 

131 ^ 


C3.57) 


iii) Qy<dy¥k: 


Wlic-n Q, < d, ¥ k, the optimal solution turns out to be 
K K 


the foil owing J 

H ' => , S"= U, = 1 ¥^ i, xj'. = 1 and w';J 


10 


1 . . ¥ J . . e S' 
iO ^0 


(3.38) 


For the production stages for which - d^^^ , the optimal 


m,' chining speeds vA^ are 




^ijk = ^ijk 


¥ J . ■ £ S 

1 J 


However, for the production stages for which Qj^ ^ ‘^k ® 
the optimal machining are determined as follows: 

We calculate the slack time, 


\ = '^k “ Sc 


82 


The optimal machining speeds are determined by formulating 
a nonlinear program of the problem 8^5 given by ( 3 . 26 ) and 

( 3 . 27 ). Hit and trial approach is used to obtain the optimal 
solution from (3.32) and (3.33). 

3 •» • 2 . 5 Optimi zing Algorith m? 

The following is the step by step description of the 
optimizing algorithm: 

Find for each stage from (3.19). If ^ k, 

{p' t) Step 2, otherwise, go to Step 3. 

S top 2 : The optima], solution is: 




Find max B. , N._ and denote the corresponding job by 

1 J J- J 


10 


J. . e R 
iO 

Set, Z = Bt"t N-r~r and go to Step 4. 
10 10 


Stpi-lA* Find a node from the nodes n s hi for which Z is 
maximum. Denote this node b 3 '' n and J-£^ in n by 
For the node n, if Z = Z, go to Step 85 otherwise, go to 
Step 5. 

Stop. .5p Find for each stage D^= h- Xr + prj^, 

max w. ^ is varied from zero to iv'v 


where ^ 

.L J 


10 ’ 10 


till <„ d^ ^ k. If wrr = 0, go to Step 6| otherwise, 


"k 


10 


go to Step 7 . 


Step C; Fiomovc the job from the set R of the node n. 

Set R = R ~ {Jo-}. If R = 0, set Z = Zj otherwise, 

find max ^ ..denote the corresponding job "by 

J^^eR ij io’ 

J rr and set Z = Z + Brr Ntt. 

13 10 10 

Return to Step 4. 

Stc-'p 7‘« Brafich the node n into two nodes, as given below. 

At first nodes Set n = n-i-1, N1 = N1 + {n}“{n}» 

For the newly created node n, set R — R~iJ^j}9 


®k 


°k = °k ^ij 


10 ' 

N ¥ J. . e R and Z = Z. 
ij 10 

N. 


If R = 0, set Z = Z; otherwise, find ^ max B^^ 
denote the corresponding job by J-j and set 
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At .second nodei Set n = n+l, N1 = N 1 + {n}. 
i'or t,he node n, set R. = Rj S = J^'y} 5 Dj^ = and 

‘^k ~ ^ Z = Z + Byy Nyy. If 

R = 09 set Z = 2; otherv/ise, recalculate the value of 

k • ^ J,. . e R and find max B. . N. .. Denote the 

^ J. . £ R 

, . 11 

corresponding job b 3 A and set Z = 2 N^-.. 

11 - 11 11 

Return to Step 4, 


Stcj)_^8^: 


The optimal solution is contained in node n. 
* 


Stpp_ 9 i 


Set S = S, Z'' = 2 , Q. 
go to Step 9. 

If Q. 


Dj^ ¥■ kj k = 1 and 


“k 


d ,.9 set V e S * and go to 




(c) 


Stop 13," otherwise, if Q, < d, , find v. and 

K h. IJK 

from ( 3 . 13 ) and (3*12) respectivel 3 ^. Calculate 

(c) 


Qi 


k 


2 

i=l 


V. 

N. 

1 

tk h 

P: 

ijk 

V. ¥ J. . £ S'^ 

Ilk 11 


X 


11 ■’ 


'*k -• ""k’ 

Step 13; otherwise go to Step 10, 

Step^lCi.- Select a job J. . e S^, Find initial value of 4 v from the 
following expressions 


^k 


(t) (c) 1/n 


V. 


-1 +(l/n^^j^-l) ( 


2^iik 


=) 


Go to Step 11. 
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xnianr * For the value of M'v obtained in Step 10, find v,. 

X ^J IV 

^2id P-jiv from (3.33) and (3»12), respectively^. 


(e) 


Calculate . 


Q. 


‘k 


Ng ^ 

3=1 ^i3 ^ 


If Q,. = dn_, set vf., = V. .. 


“k - ^k’ ''ijk ^ ''ijk 

Step 135 otherwise, go to Step 12. 




Step 12 ; Modify the value of appropriately depending on 
whether < d^ or > d^. If < d^, reduce the 
value of and if Qj^ > dj^, increase the value of p^. 
If k < Ns, set k = k+1 and go to Step 9? otherwise, 
stop. 


3.4, 2.6 Numerical 


Consider a two stage production system for which ten 
types of jobs are available for production. These jobs are classi- 
fied into four groups. The objective is to determine optimal 
machine loading and product mix when the machining speeds are 
decision variables r. The basic data for the example considered 
is given in Tables 3.3(a) and 3.3(b). 


3.4.2.? S_ olut ion; 

To find the optimal solution of the problem given in 
Sec. 3. 4. 2. 6, we follow the various steps of the optimizing 
algorithm given in Sec. 3. 4, 2. 5. The optimal solution is presented 
in Table 3.4. 
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Table 3.4: Optimal Solution of Numerical Example 6. 


Index 

of Accepted 

Quantity to be 
produced pcs 

Machining 

,m/min 

sj^eed 

Group 

Job 


At Sta^e 1 

At Stage 2, 


1 


8 

143.106 

156.739 

1 







2 


50 

133.109 

115.709 


1 


100 

127.29 

179.907 


3 


40 

110.596 

131.009 


1 


40 

163.117 

157.462 

4 







2 


50 

138.144 

171.259 



Total 

Profit = 

= Rs. 3,390.00 



Total 

Production Cost 

= Rs. 33,879.48 



Total 

Processing Time 

Utilized; 





At stage 1 

= 3,000.00 

mins* 




At stage 2 

= 2 , 500.00 

minso 


3.5 HEURISTIC SOLUTION PROCEDURES; 

In Sec. 3.4, solution methodologies for two variations 
of the problem,, namely, prespecified unit processing time and 
machining speeds are decisions variables, were developed. The 
objective was to maximize the total profit. These solution 
methodologies used the basic concept of B & B for providing 
the optimal solution when at least on one of the production 
stages, the total processing time required to process all the 
jobs to their lot sizes was greater than the available processing 


time on that stage. 
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It was realized on solving many multistage production 
system problems using the solution procedures given in Sec*3.^s 
that the computational effort increases with the number of 
the production stages and the number of jobs available for 
production. Therefore > an effort was made to reduce the compu- 
tational burden involved in the use of B & B procedure. This 
was achieved by reducing the size of the problem to be handled 
by B & B procedure. The reduction in the size of the problem is 
achieved by identifying a set of jobs which seem to have high 
potential for appearing in the optimal solution generated by 
the solution methodologi/ developed in Sec. 3. 4. 1.3. The heuris- 
tic procedure is described as follows. 

As soon as, the first step of the algorithms, described 
in Sec. 3.4, indicates that there is at least one production 
stage on which the available processing time is less than the 
total time recpired to process all the jobs to their lot sizes, 
contribution of each job towards the objective function, i.e., 
maximization of total profit is calculated. A job with the 
highest contribution is identified and is selected for inclu- 
sion in the list of selected jobs provided the available processing 
time is sufficient to process this job to its lot size. After 
the first job is selected, the next job with highest contri- 
bution is identified and the procedure is repeated till no job 
can be selected for inclusion in the list of selected jobs. 
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If the number of selected jobs in this list is less than 

three, the list is made empty and the B & B procedure, i.e.. 

Step 3 onwards of the algorithm given in Sec. 3.'^, are followed 

f ‘ 

as such. However, if the list of selected jobs contains more 
than two jobs, the last two jobs are removed from the list. 

The Step 3 of B & B procedure is suitably modified for each of 
the algorithms given in Sec. 3.4 to incorporate the heuristic 
procedure into the B & B procedure. The remaining steps of the 
algorithms, i.e.. Step 4 onwards remain the same. The following 
section gives the modified Step 3 for the various cases. 


3.5.1 Total Profit Maximization Wien Unit Processing Time is 
' Pre specified t 

3. 5. 1.1 Heuristic Algorithms 

The Step 3 of the original algorithm given in Sec. 3. 4. 1.4 


is modified as follows: 

Step 3.3^ : Select an initial node: n = 1, N1 = {n}. For this node. 


set R = U, S = 0, Sj^ = dy^ and 
= 1 ¥ i. Go to Step 3-2. 



Calculate , 


0 ¥ k, Z = 0 and 


Find max B. . N. . and denote the corresponding job 
J. .sR 

■by j^-. For each production stage, calculate 

13 

\ = + wj, where wjj = max 

w. . is varleij from zero to till 5 < i ¥ k. 

Go to Step 3 . 5 . 
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If and D, = 


^i3 " ij \ = ^k ^ S = S+{J^}, 

I *B. 1 and go to Step 8, Id = 0 or 

JijES ■^'-' 10 


'^lo ^ ^13 *' otherwise, 

set R ■ R . (J„ }, s = S + {Jj,}, Dj^ = = dj^-Dj 

and Xj = 0. 

Go to Step 3.2. 

A I If S = 0 and Wj-::. = o; set Z*= 0, S* = 0 and stop. 

-1- J 

If set S contains mor§ than two jobs, remove the 

last two jobs selected from the list S and include 

them in R ; otherwise, make S = 0, Set Xj^ = 1 for 

those groups G^, of which no job is included in S. 

Set Z = Z B ^ 1 . D = { £ (s + w p ) 

.. J. .eS ^ J. .eS 

N„ 10 10 




J, .eS 


+ X^ } and ^ Recalculate 

the value of N . . ¥ J. . e R. Find max B . . N. . and 
10 10 j_ 10 10 

denote the corresponding job by Jt-r. Set Z = Z + Bi-. N-r-» 

1(3 ^3 

and go to Step 4. 


3 . 5.2 Total Profit Maximi: 
Decision Variables s 


when Machinim 


)eeds are 


3. 5 . 2.1 


Step 3 remains the same as given in Sec, 3, 5. 1.1 

^ t ^ 

except the fact that instead of f’ijk used for 

^ r«% “II /-."I M T\T/^4- 4-V. o 4“ V / a r» c:iV^4* c* 4“ Vi ci mn o*v4 mi im, 


calculating Note that ^ 

production-rate unit processing time. 


represents the maximum- 
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3»5.3 Illustrative Example and Computational 


3. 5. 3.1 


Numerical example 5 presented in Sec. 3. 4. 1.5 is 
considered to demonstrate the heuristic solution procedure. 


3. 5. 3. 2 Solution : 

The following is the stepwise solution of the example 
using the heuristic approach. 

Step 1: = 3105.00 mins, and = 2995.00 mins. 

Since \ ^ ^ go to Step 3. 

Step 3.1s Select an initial node: set n = 1» N1 = {n} = {l}. 

For node 1, set R = '^12®' '^21' '^22” '^23®* '^31’ ^^32^ 

■^42' '^43>’ ® °k = ° ^ * >'> 

Z = 0 and = 1 V i. 

Go to Step 3.2. 

Step 3.2*. Calculate ^ '^ij ^ 

JijCR 

value is Rs. 4080.00 corresponding to job '^12* 

¥e find that = 68.00 mins and = 94,00 mins 

for w^2 = 

Go to Step 3.3* 

Step 3.3s Since ^ ^ ^^12 ~ ^12’ ^ ~ ^"’^'^12^ 

" ('^ 11 » '^ 21 * *^ 22 ’ ‘^ 23 ’ '^ 31 * ‘^ 32 * '^ 41 ’ *^ 42 * ‘^ 43 ^’ 

S = S + 

= 94.00 mins, s^ = 2332.00 mins, and s^ = 2306.00 mins. 
Go to Step 3.2. 
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Step 3. 


Now we recalculate the value of N. . ¥ J. . e R. It 

13 13 

is ohserved that max B. . N. . occurs corresponding 

to a Job J 22 * find that = 562.00 mins, 

= 473.00 mins and = 50 pcs. Go to Step 3.3. 

Since and Dj^ < d^^ ¥ k, set 

R = ‘^21’ '^23* '^31' '^32' '^41’ '^42* ‘^43^’ 

S = {J^2» '^22^’ ^2 " " 562.00 mins, = 473.00 

mins, = 1838.00 mins, and S2 = 1927.00 mins. 


The Step 3.2 and 3.3 are followed iteratively and 
sets R & S are generated, 

R = ^'^31’ '^4l* *^42^ and 

S = {J-,., J 


11^ "^IB* '^21* *^22’ *^ 23 * *^32* "^43^* Furthers 


, 3*^*1 9 3*^ 


= 2394.00 mins, D2 = 2386.00 mins, = 0 ¥ i and 


w^2 = 0- 


Since w 


42 


0 and < d^, we go to Step 3.4. 


"k 


is ¥e remove the last two jobs J2J and from S 
and include them in R. 


Now 5 set R — ^'^23^ '^31’ ^^41* "^42’ *^43^ ^ 


(J- 


J-IO. J , JoO» 3''2,p} , X 


■4 ~ i» 


'll’ "12’ "21’ "22* 32 

Z = Rs, 2210.00, = 1504.00 mins, 

D2 = 1536.00 mins, s^ = 896.00 mins and 

S2 = 864.00 mins. 
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Recalculate N. . ^ J. • e R. Find max B. . N. .. Its 

la la 11 11 

value is Rs. 828,00 correspcnding^'fio job ^ 23 * 

Set 2 = Z + B 2 ^ N 2 ^ = Rs. 3038,00 and go to Step 4, 

Now we follow B & B procedure from Step 4 onwards as 
discussed in Sec, 3. 4, 1,4, The final product mix obtained is 
as shown in Table 12. 

3 • 5 . 3 . 3 C omput a tional experience s 

It was observed that the solution of the numerical 
example 5 obtained b 3 r heuristic and optimizing procedures turns 
out to be the same. The heuristic and optimizing solution 
procedures required a CPU time of 0,20 secs, and 0.26 secs, 
respectiveli/' for solving this example ■■■nDEC-1090 computer 
s'/stem. Number of problems were solved using both the approaches. 
In each case, it was observed that the heuristic approach requires 
less computational effort as compared to the optimizing procedure. 

3.6 VARIATIONS OF MULTI-STAGE PROBLEM J 

In Sec, 3.4 and 3.5, we presented mathematical 

models and solution procedures for the multi-stage production 

systems and the objective was to maximize the total profit. The 

mathematical models and solution procedures can be suitably 

modified for maximization of total number of units produced. 

The term B. i.e,, the per unit profit for job J. is 
1 J J 

replaced by 1. When Ns is put equal to 1 , the multi-stage 
problem becomes single stage problem. 
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The multi-stage algorithms given in Secs. 3.^ and 
3.5 were used to solve the single stage numerical examples 
1,25 3 and 4 and it was observed that the solutions obtained 
were same as obtained by algorithms given in Sec. 2.4. 

The solution procedure given in Sec, 3.4.2 is based 
on the fact that the unit processing times of the jobs are 
determined by machining speeds which are treated as decision 
variables. However, there may be multi-stage problems in 
which the unit processing time for some jobs are prespecified 
while for the remaining jobs the unit processing time depend 
upon the machining speeds which are also to be optimally 
selected. 

Initially, we assume that the jobs for which the 
machining speeds are not prespecified, are processed using 
minimum unit production time machining speeds. The machine 
loading and product mix is determined using the algorithm 
given in Sec. 3.4.1. Using this product mix, the time required 
on each production stage to process the selected jobs for 
which the unit processing time is not prespecified, is calcu- 
lated and subtracted from the available processing time for that 
stage. The remaining available processing time at each stage 
is used to obtain the optimal machining speeds for the various 
jobs (for which the processing times are not prespecified) to 
be processed on that stage. 
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The multi-stage problem is now converted into 
multiple single stage problems (the number of single stage 
problems considered being equal to the number of stages in 
the multi-stage problem) . The optimal machining speeds are 
determined using the basic concepts given in Sec. 2. 4. 1.2. 



CHAPTER IV 


CONCLUSIONS AND SCOPE FOR FURTHER WORK 

4.1 CONCLUSIONS; 

In this thesis, we have developed mathematical models 
and solution procedures for the machine loading and product 
mix decision problems for single and multi-stage production 
systems with limited available production times at each stage. 
The objective was to maximize either the total profit or the 
total amount of production. Two important variations of the 
problem, namely, the unit production time is prespecified and 
machining speeds are decision variables, were also considered. 

For each one of the above stated variation of the prob- 
lem, an optimizing algorithm which utilized the basic frame- 
work of branch and bound procedure, was developed. In order to 
reduce the computational effort, heuristic algorithms were 
also developed for each of the problem variations. Sample 
problems were solved using both the optimizing and heuristic 
algorithms. 

4.2 SCOPE FOR FURTHER WORK; 

In order to account for realistic production situations, 
the single and multi-stage models developed in this thesis need 
to be extended to account for the following situations; . 
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i) Each production stage has one or more machines. 

ii) Machines are capable for processing more than one unit 
at a time. 

iii) The jobs may have technological ordering for the opera- 
tions to be performed on various production stages. 

iv) The jobs may have due dates associated with them. 

In Section 2.4. 1.2, a technique was suggested for deter- 
mining the step length for decreasing or increasing the value 
of jj, . Better methodology^ needs to be developed for appro- 
priate selection of step lengthj since this will considerably 
reduce the computational effort involved in the use of various 
algorithms. 
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